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PREFACIO 


En febrero de 1979 se firmo un convenio de colaboración entre la 
UNAM, PEMEX, IMP y el CIPM (Colegio de Ingenieros Petroleros de 
México). El objeto del convenio ha sido elevar el nivel academi. 
co de los alumnos del área de Ingeniería Petrolera en la Facul — 
tad de Ingeniería, tanto de licenciatura como de posgrado, así - 
como crear el Doctorado, y promover la superación de un mayor nú 
mero de profesionales que laboran en la industria petrolera, por 
medio de cursos de actualización y especiali2ación . 

Uno de los programas que se están llevando a cabo a nivel de li- 
cenciatura, dentro del marco del Convenio, es la elaboración y - 
actualización de apuntes de las materias de la carrera de Inge — 
niero Petrolero. Con esto se pretende dotar al alumno de más y - 
mejores medios para elevar su nivel académico, a la vez que pro- 
porcionar al profesor material didáctico que lo auxilie en el — 
proceso enseñanza-aprendizaje. 


Estos apuntes fueron preparados utilizando en parte material con 
tenido en el libro "FLUJO MULTIFASICO EN TUBERIAS", elaborado — 
por los ingenieros Antonio Acuña R. , Francisco Garaicochea P. y 
Tomás Limón H. ( y publicado por el IMP en 1976. 

Con el objeto de facilitar el empleo de computadoras programa- 
bles, se realizó una investigación sustancial para presentar — 
ecuaciones, en vez de las gráficas que muestran los resultados 
de las correlaciones. 


Es conveniente indicar que los temas tratados en estas notas — 
complementan a los elaborados por el ingeniero Eduardo Lozano - 
V. , en sus "APUNTES DE MEDICION Y TRANSPORTE DE HIDROCARBUROS". 


Se agradece la colaboración del estudiante Miguel Angel Hernán- 
dez García, en la revisión de algunos capítulos de estos apun — 
tes. 


DEPARTAMENTO DE EXPLOTACION DEL PETROLEO 


Noviembre de 1983. 
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1.- PROPIEDADES DE LOS PUJIDOS 


1.1. - Introducción. 

El primer problema que surge en relación con la determinación de las pro 
piedades de los fluidos, es la carencia de análisis pVT apropiados de la 
boratorio. El análisis con que generalmente se cuenta está realizado a “ 
la temperatura del yacimiento, bajo condiciones de separación diferen- 
cial^sin embargo, al escurrir los fluidos a través de la tubería de pro 
ducción y de la línea de descarga, su temperatura disminuye y el gas li~ 
berado no se separa totalmente del aceite. Para conocer, a diferentes pre 
siones y temperaturas, las propiedades de los fluidos, estas se determi-~ 
nan generalmente por medio de correlaciones. Al usar correlaciones se so 
breentiende que se obtendrán valores aproximados de las propiedades men- 
cionadas, lo que en sí acarrea un margen de error. 

Para facilitar el uso de calculadoras programables, los resultados de las 
correlaciones se expresan en forma de ecuaciones, en lugar.de presentar 
las figuras que generalmente aparecen en los trabajos originales. 

Las unidades en que están expresadas las ecuaciones son las que aparecen 
en la nomenclatura, a menos que se indiquen en el texto otras unidades. 

1.2. - Propiedades del aceite saturado. 

^ 1 * 

1.2.1.- Correlación de M.B. Standing. 

Esta correlación establece las relaciones empíricas observadas entre la 
presión de saturación y el factor de volumen del aceite, en función de 
la razón gas disuelto - aceite, las densidades del gas y del aceite pro 
ducidos, la presión y la temperatura. La correlación se estableció para 
aceites y gases producidos en California y para otros sistemas de crudo 
de bajo encogimiento, simulando una separación instantánea en dos eta- 
pas a 100 °F. La primera etapa se realizó a una presión de 250 a 450 - 
lb/pg 2 abs, y la segunda etapa a la presión atmosférica. 

Debe entenderse que la densidad del aceite producido en el tanque de al- 
macenamiento dependerá de las condiciones de separación (etapas, presio 
nes y temperaturas) . Mientras más etapas de separación sean, el aceite - * 
será más ligero (mayor densidad API) . 

La presión del aceite saturado se correlacionó en la siguiente forma: 


R 0.83 

( — !— > 

10 0 .00091(T) "| 

F g 

iqO-o^Oo’J 


Por lo que despejando la relación gas disuelto - aceite (R ) de la ecua 
ción anterior se tiene: 


* Referencias al final de cada capítulo. 


Vn. opZtdadu dt lo¿ Fluidos . 


R 

s 


El factor de volumen del aceite fue correlacionado con la relación gas 
disuelto-aceite, la temperatura, la densidad relativa del gas y la den- 
sidad del aceite. Se obtuvo la siguiente expresión: 

1 175 

B = 0.972+0.000 147 (F) (1.2.3) 

o 

Donde : • 

F = R ( y / y ) 1/2 + 1.25 T (1.2.4) 

s r g r o 


. 2 

1.2*2.- Correlación, de Vázquez 

Para establecer estas correlaciones se usaron más de 6000 datos de R^, 

B y p a varias presiones y temperaturas. Como el valor de la densi- 
dad relativa del gas es un parámetro de correlación importante, se deci- 
dió usar un valor de dicha 2 densidad relativa normalizado a una presión 
de separación de 100 lb/pg manometrica . Por lo tanto, el primer paso 
para usar estas correlaciones consiste en obtener el valor de la densi- 
dad relativa del gas a dicha presión. Para esto se propone la siguiente 
ecuación : 

v = v i 1 + 5.912 x 10**" 5 y T log (p / 114.7)). (1.2.5) 

p gs r gp -o s ^ s 


densidad relativa del gas resultante de una separación a 100 
lb/pg 2 manometrica. 

densidad relativa del gas obtenida a las condiciones de sepa- 
ración de p g y T g . 

- 2 

presión de separación real, en lb/pg abs 
temperatura de separación real, en *F. 

La correlación para determinar R se afinó dividiendo los datos en dos — 
grupos, de acuerdo con la densidad del aceite. Se obtuvo la siguiente - 
ecuación: 

R = C v p C 2 exp (r ( y A T+ 460))) (1.2.6) 

s 1 r gs 3 o 


Donde : 



T 

s 


[ 


18 


10 0.0125 ( r o ) 

0.83 

10 0 .0009.1 (T) 



( 1 . 2 . 2 ) 
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áT ^5^5*^ 


Los valores de los coeficientes son: 


Coeficientes 

C, 


Y Q — 30 ° API 


0.0362 

1.0937 

25.724 


r >30° API 
o 


0.0178 

1.1870 

23.931 


La expresión que se obtuvo para determinar el factcr de volumen es: 

B o =1 + C 1 R s + c 2 (T ' 60) <W + c 3 R s (T - 60) ( r o / V gs ) (1.2.7) 


Los valores de los coeficientes son: 


Coeficientes 


X ^ 30 ° API 
o 


4.677 x 10 


-4 


-5 


1.751 x 10 


-1.811 x 10 


V q >30° API 


4.67 x 10 


-4 


-5 


1.1 x 10 


1.337 x 10 


1.2.3.- Correlación de 0istein 

Esta correlación fue establecida, utilizando muestras de aceite produci- 
do en el Mar del Norte, donde predominan los aceites de tipo volátil. 
Los valores de y B q se obtienen mediante los pasos siguientes : 

* 

1». - Calcule p con: 

log p* = - 2.57364 + 2.35772 log p - 0.703988 (log p) 2 + 0.098479 

(log p) 3 . . (1.2.8) 


2 . - Calcule R con: 


R * y 

s l g 


- ' 0,989- 

' p V ' 

) 

m 0.130 1 

L T J 


i 


0 .*816 


3 . - Calcule B * con : 


0.526 


B * = R ( X / X ) 

o s g ro 


+ 0.968 T 


(1.2,9) 


( 1 . 2 . 10 ) 
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4.- Determine B con: 
o 


log (B - 1) = - 6.58511 + 2.91329 log B - 0 .27683 (log B ) (1.2.11) 


1.2.4.- Densidad del aceite saturado. 


La densidad del aceite saturado, en Ib /pie 3 , es: 

m 


* 62.4 0.01362 Rs Y ad 


( 1 . 2 . 12 ) 


1.2.5.- Correlación para determinar la viscosidad del aceite saturado. 4 

La viscosidad del aceite saturado se puede calcular de la manera siguien- 
te: 


"o = a 

(1.2.13) 

a = 10.715 (R + 100) - °' 515 
s 

(1.2/14) 

b = 5.44 (R + 150) "°‘ 338 
s 

(1.2.15) 

M om - 10* - i 

(1.2.16) 

X = y T - 1 - 163 

(1.2.17) 

y =■ 10 Z 

(1.2.18) 

Z = 3.0324 - 0.02023 Y 
o 

(1.2.19) 


Donde: 


: Viscosidad del aceite saturado 


^om 


: Viscosidad del aceite muerto a T 


1.2.6.- Correlación para calcular la tensión superficial del aceite 
saturado . ^ 

La tensión superficial del aceite saturado, en dinas /cm, se puede deter- 
minar con la siguiente expresión: 


a Q = (42.4 - 0.047 T - 0.267 V\) exp (-0.0007 p) 


( 1 . 2 . 20 ) 
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1.3.- Propiedades del aceite baj osaturado. 

1.3.1.- Correlación para obtener la compresibilidad del aceite baj osaturado/ 

La ecuación siguiente sirve para determinar la compresibilidad del acei- 
te baj osaturado : 


(a l + *2 R s + a 3 T + a 4 ''gs + a 5 V 7 a 6 P 


(1 . 3.1) 
5 


Donde: a 1 = - 1433, = 5, a 3 = 17.2, a 4 = - 1180, a 5 = 12.61, a & = 10 

1.3.2.- Densidad del aceite baj os aturado 

La densidad del aceite bajosaturado está dada por la siguiente errpre — - 
sión: 


Po “ p ob exp C ob ( P - P b > 


(1.3.2) 


Donde : 


: densidad del aceite a 
C ob 1 compresibilidad del aceite a p fc 

fe ^ 2 

1.3. 3. - Correlación para calcular la viscosidad del aceite bajosaturado. 
La viscosidad del aceite bajosaturado se obtiene de la manera siguiente: 


^ob 


<p/p h ) 


= C T p 2 exp (C '+ C p) 


(1.3.3) 

(1.3.4) 


Donde: 


= 2.6 


C 2 = 1.187 


11.513 


C . * - 8.98 X 10 
4 


-5 


"ob ‘ viscosidad a P b 

1.3.4.- Correlación para obtener el factor de volumen del aceite 
2 

bajo saturado. 


Paira el aceite bajosaturado se tiene la ecuación: 
B B , / exp (C (p - P, ) ) • * * * 

O — OD O XD 


(1.3.5) 


Donde Co es la compresibilidad del aceite bajosaturado, calculable con 
la Ec. (1.3.1). 
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1.3.5.- Correlaciones para obtener p, 

b 

Es obvio que el primer paso para obtener las propiedades del aceite bajo- 
saturado es la determinación de la presión de saturación del aceite. 

Esta presión puede obtenerse con la Ec. (1.2.1) cuendo el aceite es de - 
bajo encogimiento. Para aceites volátiles se puede usar la correlación 
de 0istein: 


Donde : 


log p b - 1.7669 + 1.7447 log p * - 0.30218(log p*) 2 (1.3.6) 


* / R x 


0.816 


m 0.130 
T 

y 0.989 


(1.3.7) 


El valor de p^ así obtenido puede corregirse usando una ecuación similár 
a la (1.19). Es conveniente indicar además que generalmente la presión 
de saturación del agua se considera igual a la del aceite. 


(12) 

1.4.- Propiedades dal gas natural 
1.4.1.- Densidad relativa del gas 

Generalmente se utiliza sólo el valor de la densidad relativa del gas 
producido (que se proporciona como dato) en los cálculos de las propie- 
dades de los fluidos. Sin embargo, es conveniente distinguir la densi- 
dad relativa del gas libre ( y ) y la densidad relativa del gas disuel- 

gf 

to ( y . ) de la densidad relativa del gas producido ( y ) . 
gd g 

La densidad relativa del gas disuelto puede obtenerse con: 

y ■ = 0.25 + 0.02 y + (10 ~ 6 ) (0.6874-3.5864 y ) R . . . . (1.4.1) 

gd o o s 

El valor de la densidad relativa del gas libre se obtiene con: 

y ~ = (R Y -R y ,) / (R - R ) . . . . . . . . - (1.4.2) 

r gf f g s r gá ' s 
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i. 4.2.- Determinación del factor de volumen del gas. 
De la ecuación de los gases reales se obtiene; 


B = Q. 02825 Z (T + 460) 
g P 


(1.4.3) 


1.4.3.- Determinación de la densidad del gas. 

La densidad del gas esta dada por la siguiente expresión; 

p = 0.0764 Y af 

g B 


(1.4.4) 


Substituyendo (1.4.3) en (1.4.4) 


= 


2.7044 P 


Z (T + 460) 


(1.4.5) 


1.4.4.- Determinación del factor de compresibilidad del gas/'^ 

Existen diferentes correlaciones para calcular las propiedades fy pseudocrí- 
ticas del gas. Las curvas correspondientes a los gases han sido estable- 
cidas utilizando gases de los separadores y vapores obtenidos en los tan- 
ques de almacenamiento. Estos gases contienen altas proporciones de meta- 
no, etano, propano y butano. Las curvas correspondientes a los "condensa- 
dos" pertenecen a gases que contienen cantidades relativamente grandes de 
los componentes mas pesados. Standing^, sugiere el uso de las curvas de 
"condensados" en los cálculos que involucren gases en equilibrio con el - 
aceite y el uso de las curvas correspondientes al gas , para gases super- 
ficiales . 

La ecuación para gases superficieales es ; 


T = 167 + 316.67 y 
pe ' gf 


p = 702.5 
PC 


50 V 


gf 


(1.4.6) 

(1.4.7) 


La ecuación para gases húmedos (condensados) es: 


pe 


238 + 210 T 


gf 


(1.4.8) 


p = 740 - 100 V 
pe gf 


(1.4.9) 


Las ecuaciones siguientes permiten calcular, por ensaye y error, el valor 
de Z, usando como datos: 


T = 
pr 


T + 460 


pe 


(1.4.10) 
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P r 


pe 


(1.4,11) 


p = 0.27 P / Z T 
r pr pr 


(1.4.12) 


Z = (A 1 + A 2 / T pr + A 3 / T pr ) P r > CA 4 + A 5 /T pr ) ^ 


+ (A 5 A 6 p r 5 ) / T Pr + (A, P r 2 /T pr 3 ) (1 + AgP/) 

(exp (-A P 2 )) + 1 
o r 


(1.4.13) 


Donde : 


A ± * 0.31506 


A 2 


= - 1.0467 


A 3 = - 0.5783 
A, - 0.5353 


A 5 ” 


A 7 * 


- 0.6123 

- 0.10489 
0.68157 
0.68446 


El procedimiento consiste en suponer un valor de Z y obtener p » para 
ese valor supuesto. Se calcula Z con la Ec. (1.4.13) y se compara con 
el supuesto. Si no coinciden estos valores, se supone para la s iguien 
te iteración el valor de Z calculado. El procedimiento se repite has_ 
ta caer dentro de una tolerancia preestablecida (menor : o igual a 0.001) 


1.4.5.- Propiedades de gases que contienen N 2 , CC> 2 y H 2 S 
Las propiedades pseudocríticas de gases que contienen cantidades apre- 
ciables de N 2 , C0 2 y H 2 S, pueden calcularse por el método de Wichert y 
Aziz. 10 La densidad relativa de este gas puede calcularse con: 

V gc = ( T g ” 0*967 yN 2 - 1.52 yC0 2 - 1.18 y^S) 

/a- yN 2 - yC0 2 - yH 2 S) U.4.14) 


las propiedades pseudocríticas de dicha mezcla se obtienen con: 


pem 


= (1 - yN_ - yCO - yH S) P +493 yN, 


+ 1071 yC0 2 + 1306 yH 2 S 


pe 


(1,4.151 


T ^ d “ YN, 

pem ¿ 


yC0 2 - yH 2 S) 


T + 227 yN„ 
pe 1 2 


+ 548 yC0 2 + 672 yH 2 S 


(1.4.16) 
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Estos valores deben de ajustarse, por su contenido de co.'y H S, median 
te las ecuaciones de Wichert y Aziz: 2 2 

C wa = 120(( y C0 2 +YH 2 S)0 ’ 9 - ( y co 2 + y H^S) 1 - 6 ) 


+ 15 (y H 2 S°- 5 - y ¡^S 4 ) 


T' * T - C 
pcm pcm wa 


(1.4. 17) 

(1.4.18) 


P 'pcm = P pcm T 'pcm 7 (T pcm + y H » síl " y H , s > C .,.J (1-4.19) 


2 wa 


1.4,6.- Correlación para determinar la viscosidad del gas. 

La Viscosidad del gas se obtiene con la siguiente ecuación: 7 


/u g = K (*0”Vexp£(X) (Ai/62.428 ) y ] ' U.4.20) 


(9.4 + 0.5794 V Qf )( T + 460) 


1.5 


209 + 550.4 y + (T + 460) 
gf 


(1.4.21) 


X - 3.5 +986 


_ + 0.2897 y 

(T + 460) gf 


(1.4.22) 


Y = 2.4 - 0.2 X 


(1.4.23) 


1.5.- Propiedades del agua saturada: 

1.5.1.- Correlación para determinar el factor de volumen del agua satu- 
rada. 8 

El factor de volumen del agua saturada se puede calcular con la 
siguiente ecuación: 

B = 1.0 + 1.2 X 10” 4 (T-60) + 1.0 x 10 " 6 (T-60) - 3.33 x 10~ 6 p 


(1.5.1) 


1.5.2.- Densidad del agua saturada. 


La densidad del agua saturada puede obtenerse de la expresión 
siguiente: 

62.43 


^w= B 


(1.5.2) 
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1.5.3. - Correlación para determinar la viscosidad del agua saturada. 9 

La viscosidad del agua saturada es funci6n de el porcentaje de NaCl que 
contenga y está dada por: 

% - A + B/T (1.5.3) 

Donde: 

A - - 0.04518+0.009313 (% NaCl) - 0.000393 (% Na Cl) 2 (1.5.4) 

B= 70.634 + 0.09576 (% NaCl) 2 (1.5.5) 

Si las presiones son elevadas, es necesario corregir el valor de la vis 
cosidad, obtenido con la Ec. (1.5.3), por efecto de la presión. Este 
factor se obtiene con la expresión: 

f(pr T) * 1 +3.5 x 10“ 12 (p) 2 (T-40) (1.5.6) 

Esta correlación puede aplicarse para el rango de valores siguiente: 

60 °F T ^400°F 
P ^ 10 000 lb/pg 2 
Salinidad % Nací ^ 26 % 

1.5.4. - Correlación para calcular la tensión superficial agua - gas. 

La tensión superficial agua - gas se calcula con las siguientes expre- 
siones: 


*wl ■ 

52.5 - 0.006 p 


(1.5.7) 

°w2 '* 

76 exp (- 0.00025 p) 


(1.5.8) 


1 (280 - T ) 

1 


* 

206 °w2 ~ °wl \ + 

<*wl 

(1.5.9) 

^wl : 

: tensión superficial agua - gas 

a 280 °F 


°w2 : 

¡ tensión superficial agua - gas 

a 74 °P 


°w ! 

:* tensión superficial a p y T 
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1.5.5.- Correlación para determinar la solubilidad del gas en el agua. 


La R se calcula de la siguiente manera: 
sw 

i 

p = 1 - exp (- p/2276) 

T* = (5/9) (T - 32) 

_ T* - 90 
10 

S - P ' (A + BT ' + CT ' 2 + DT ' 3 ) 

Donde: 


A = 3.69051 
C * 0.01129 


B = 0.08746 
D = -0.00647 


R SW - 5 . 6146 S . 


(1.5.10) 

(1.5.11) 

(1.5.12) 

(1.5.13) 

(1.5.14) 


El valor de R¿w así obtenido, debe corregirse para considerar el 
efecto de la salinidad del agua. El factor de corrección es 

f cs 3 1 + (0.0001736 - 0.Q7703 % NaCl ) (1.5.15) 

1.6,- Propiedades del agua bajosaturada 

8 

1.6.1.- Compresibilidad del agua bajosaturada . 

La compresibilidad del agua se puede determinar de la siguiente manera: 

( 1 . 6 . 1 ) 


C = (A + BT + CT 2 ) 10 " 6 f* 
w 


donde : 


A =3.8546 -0.000134 p 


*1 


B = - 0.01052 + 4.77 x 10 
C = 3.9267 X lo" 5 - 8.8 x lo" 10 p 


f* = 1 + 8.9 x 10 


-3 


( 1 . 6 . 2 ) 

(1.6.3) 

(1.6.4) 

(1.6.5) 


Donde : 


f* = factor de corrección por presencia de gas en solución. 


i-ir 
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El rango de aplicación de esta correlación es: 


1000 Ib/pg abs ¿lp 


6000 lb/pg abs 


80 


. 250 °F 


3 3 

0 pies /bl >1 R gw ^ 25 pies /bl 


1.6.2.- Factor de volumen del agua bajosaturada. 

El factor de volumen del agua bajosaturada esta dado por la siguiente 
expresión: 


B 


) 


-W B wb exp ( - C «b (p - V 

1.6.3.- Densidad del agua bajo saturada. 

La densidad del agua bajosaturada se determina con la ecuación 


( 1 . 6 . 6 ) 


P v m P vb e *P ( - c „b ( P - P b > 1 (1.6.7) 

1.7.- Ejemplos de aplicación: 

Ejemplo 1.1. 

Un pozo de exploración produce aceite de 30°API con un relación gas- 
aceite instantánea de 358 pies 3 /bl. Se mide su presión estática que es 
de 3400 lb/pg 2 abs. La densidad relativa del gas producido es de 0.75. 
El aceite a condiciones de yacimiento ¿Es taí saturado o baj osaturado? 
si la temperatura del yacimiento es de 200 °F. 

Solución : 


Se determinará la presión del aceite saturado, considerando el valor de 

!. ( 1 . 2 . 1 ) 

0,00091 (T) 


R como el de R g . Aplicando la Ec. (1.2.1) 


p * 18 ( 


p = 18 1 

Para este caso p = p^ 


[■ 

‘ 1 0.75 1 


10 


10 


0.0125 ( V c 


10 


0.00091(200) 


10 


0.0125(30) 


-] 


por lo que se obtiene : 


p, » 1930.61 lb/pg abs 
b 

De el resultado obtenido se concluye que el yacimiento es bajos aturado 
ya que . < P ws 


1-12 


Propiedades de. lo * fluido* 


Ejemplo 1.2.- 

E1 siguiente problema se resolverá usando las tres correlaciones que se 
vieron para determinar las propiedades del aceite sat arado. El ejemplo 
corresponde al de un aceite volátil producido por Pemex en la Zona Sur. 
Los resultados se compararán con los reportados de i laboratorio para ver 
cual de las correlaciones utilizadas da valores mas aproximados. 

Debe indicarse que el análisis pVT obtenido en el laboratorio general- 
mente corresponde al de un proceso de separación diferencial, mientras 
que las correlaciones que aquí se ven corresponden a un proceso de sepa- 
ración instantánea. 

Se desea determinar las propiedades de un aceite saturado a una presión 
de 4500 lb/pg 2 , abs, a una temperatura de 240 °F si V = 38°API, . Y =0.8. 
El aceite producido, de acuerdo a sus características? se cons idera^como 
probablemente de tipo volátil. Las condiciones de separación fueron.. 

P s =120 lb/pg 2 abs y T s = 80 °F. 

Utilizando la coorrelación de Standing Ec. (1.2.2) 



R 


s 


1263 



Utilizando las Ecs. (1.2.4) y (1.2.3) se calcula el factor de volumen 


F = R ( V / V) 
s g o 


F= 1263 (0.8/38) 1/2 + 1.25 (240) = 483.255 
1 175 

B = 0.972 + 0.000147 (F) 
o 

B = 0.972 + 0.000147 (483.255) 1 - 175 


B = 1.1815 

o 
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Se observa que el factor de volumen del aceite obtenido, al aplicar la 
correlación de Standing, resulta muy bajo y no corresponde al de un - 
aceite supuestamente volátil. Esto se debe a que la correlación de Stan- 
ding fuá desarrollada utilizando muestras de aceite de bajo encogimien- 
to. Por lo tanto, estos resultados (los de R y B ) no son confiables. 
Será necesario utilizar otras coorrelaciones S para°calcular los valores 

de R y B . 
s o 

Aplicando la correlación de Vázquez se calcula la densidad del gas nor- 
malizada a una presión de 1001b /pg 2 con la Ec. (1.2.5). 

10“ 5 y T log (P /114.7) 
os s 

10~ 5 ( 38} 80 log £120/114.7) 

la relaciSn de solubilidad (R s ) 

Y /(T + 460))) 

O 

1 187 

R g = 0.0178 (0.8028) 4500 exp (23.931(38/(240 + 460)) ) 

R = 1136.54 pies 3 gd a c,s, 

bl a c.s. 
o 

La determinación del factor de volumen se realiza con la Ec. (1.2.7) 

B = 1 + C. R + C_ (T-60) ( yj y ) + C-R (T-60) ( }’ / y ) 
o 1 s 2 o gs 3 s o f gs 

B = 1 + 4.67 x 10 -4 (1136.6) + 1.1 x lo” 5 (240-60) (38/0.8028) 
o 

+ 1.337 X 10“ 9 (1136.6) (240 - 60) (38/0.8028) 

B a 1.6245 
o 

Con la correlación de 0i»tein se tienen los siguientes resultados: 

De la Ec. (1.2.8) 

2 3 

log p* =. - 2.57364 + 2.35772 log p - 0.703988(log p) + 0.098479 (log p) 

log p* = - 2.57364 + 2.35772 log (450p) - 0.703988 (log 4500) 2 
+ 0.098479 (log 4500) 3 


y « y (1 + 5.912 x 
gs gp 

y = 0.8 (1 + 5.912 X 
gs 

v = 0^028 

r gs 

Con la Ec. (1.2.6) se calcula 

R = C y p C2 exp (C ( 
s 1 gs 3 
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p* = 27.901 

Con la Ec. (1.2.9) se obtiene R 


R = V 


0.989 


0 . 130 


0.816 


R = 0.8 
s 


27.901 (38) 


0.989 


(240) 


0.130 


0.816 


R = 1622.54 pies 3 g.d. a c.s. 

s bl a c.s. 

o 


La densidad relativa del aceite está dada por: 

= 141.5 

y ro 131.5+'° API 

Por lo que para el ejemplo. 


141.5 


F ro “ 131.5 + 38 


= 0.835 


Con la Ec * (1.2.10) 


B = R ( 
o s 


0.526 

) + 0.968 T 


V = 1622,-5412 <-—3 5 


0.526 

) +0.968 (240) 


Bo = 1819 

Se determina el factor de volumen con la Ec* (1.2.11) 


* . 2 


log (B - 1) = - 6.58511 + 2.91329 log B q - 0.27683 (log B 0 ) 

log (B - i)= - 6.58511 + 2.91329 log 1819 - 0.27683 (log 1819)" 

log (B q - 1) = - 0.0300085 

B = 1.9332 
o 
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A continuación se muestra una tabla comparativa de los valores obtenidos 
con cada una de las correlaciones utilizadas y el reportado de laborato- 
rio . 


Correlación de Standing 
Correlación de Vázquez 
Correlación de 0istein 
Laboratorio 


R 

B 

s 

o 

1263 

1.1815 

1137 

1.6245 

1623 

1.9332 

1725 

1.9402 


Los valores calculados con las correlaciones pueden corregirse para que 
concuerden con los obtenidos en el laboratorio. Para ésto se usa el pro- 
cedimiento siguiente, aplicado a R : 


R (R . . /R ) 
se slab se 


(1.19) 


Donde : 


R : valor de R corregido 

sr s 

R : valor de R obtenido con la correlación 

se s 

r : valor de R obtenido del laboratorio, 

slab s 


(R /R ) : factor de corrección, 

slab se 


Ejemplo 1.3.- 

Calcular la densidad y la viscosidad del aceite del ejemplo anterior 
utilizando los valores obtenidos en el laboratorio. 

Con la Ec. (1.2. 12) se obtiene el valor de la densidad: 

62.4 7 + 0.01362 R V 

p ro s q 

° = 

_ 62.4 (0.835) + 0.01362 (1725) (0.8) 

o " 1.9402 

p - 36.54 Ib /pie 3 
o 
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Con las Ecs. ti. 2. 13) a la (1.2.19) se calcula el valor de la viscosi- 
dad. 

Z =3.0324 - 0.02023 V 
o 

Z = 3.0324 - 0.02023 (38) = 2.2636 


.-Z , _2.2636 _ __ 

y = 10 = 10 = 183.4847 


x = y T 


-1.163 


X = (183.4847) (240)” 1 ' 163 = 0.3129 


Pcm 

li 

s 

- 1 


i -0 , 

. 3129 

Posa 

= 10 


que es la viscosidad del aceite de densidad 38°API, sin gas disuelto, 
a una; temperatura de 240 °F . ‘ 


b = 5.44 (R + 150) 
s 


-0.338 


b = 5.44 (1725 + 150) ’ 338 = 0.4259 


a = 10.715 (R + 100) 
s 


- 0.515 


a = 10.715 (1725 + 100)“°* 515 = 0.2241 


Po 


Pom 


P o 


0.2241 (1.05565) 


0.4259 


Po 


=0.2293 cp 


que es la viscosidad del aceite, de densidad 38°API , con gas disuelto, 
a Una presión de 45001b/pg^ y a 240 °F 
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Ejemplo 1.4.- 

Calcular la compresibilidad del aceite del ejemplo 1.2 a la presión de 
saturación = 4500 lb/pg 2 abs y a una presión de 6000 lb/pg 2 abs : ade- 
mas, a ésta presión, calcular la densidad y la viscosidad. Utilice los 
valores de R y B que se obtuvieron del laboratorio. 

Utilizando la Ec. (1.3.1) 


(a l + a 2 R s 


a_ T + a. 
3 4 


C = (-1433 + (5) (1725) + (17.2) (240) + (-1180) (0.80-28) + (12.61) (38) ) 

° 5 

/ 10 (4500) 


El valor de F se obtuvo en el ejemplo 1.2 con la coorrelación de 
Vázquez. 

C « 2*4115 X 10~ 5 pg 2 / Ib a p * 4500 lb/pg 2 y T * 240°F 
o 

C « (- 1433 + (5) (1725) + (17.2) (240) + (-1180) (0.8028) + (12.61) (38) ) 
o 

/ 10 5 (6000) 



Substituyendo este valor en la Ec . (1.3.2) 

P o = p ob exp (C ob <P - Pb” 

P =36.542 exp (2.4115 X 10 _S (6000 - 4500)) 
o 

p a 37.888 lb/pie 3 
o 
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2 

que es la densidad del aceite a p = 6000 lb/pg y T = 240 °F 

Utilizando las Écs. (1.3.3) y (1.3.4) se obtiene el valor de la viscosi 
dad. 

C2 

m *■ C 1 p exp CC 3 + C 4 p) 

m (2.6) (6000) 1,187 exp (-11.513 + (-8.98 X 10 “ 5 (6000)) 
m = Q .4630 

K" m . M 0b <p/p b i m 

utilizando el valor de jtí^ calculado en el ejemplo 1.3 . 
u = (0.22933) (6000/4500) 0,4630 

O 

fi Q ss 0.2619 c.p . a p = 6000 lb/pg 2 
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FIG . 2.1.- DIAGRAMA DE FLUJO EN UN CONDUCTO AISLADO. 
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Fundamento* de Flujo a Tsiavéó de TubeJita * 


De acuerdo con la ley de la conservación de la energía; 


E. + AW,. + Aw 
1 f s 


( 2 . 1 ) 


Donde : 


Á w f -“ perdidas de^energía por fricción. Estas pérdidas corresponden 
a la fricción interna del fluido (viscosidad) y a la fricción 
del fluido con las paredes rugosas de la tubería. 

4W perdidas de energía por trabajo externo. 

E^ . - energía por unidad de masa, en la posición 

,1b ¿r - pie , 
uno ( r _ ) 

Ib 


E 2 .- energía por unidad de masa, en la posición 

dos ( lb r ¿ - P le ) 

m 


Energía de expansión tÉ ) 


La energía de expansión está dada por; 


r. / ^f “ „ /l bf 

Ee { lb~ * " p 

m pie^ 


-) V (gs_) =pv. 


( 2 . 2 ) 


Donde; 

. 3 

V.- volumen específico — ) 

Ib 

m 

Energía potencia l(Ep ) 


La energía potencial está dada por la expresión siguiente : 

2 

, lbf - pie v , pie , 1 , lbf - seg , , 

E p (-jg £- ) = g ( g-) ( ,^1 — — ¿Z •) h 


2 g 
seg ^c 


Ib - pie 


(pie) 


( 2 . 3 ) 


Energía Cinética (Ec) 


La energía cinética se expresa como sigue; 

.2 ,2 


.2 


EC t i* --f (-EÍVi f < ffi - ' - S- - ’ - ir •• <’•« 

seg ^c m - * • c 


Ib 
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Donde : 


v. - velocidad 

seg 


Al sustituir las energías correspondientes a las posiciones 1 y 2 en la 
ecuación (2,1) se obtiene: " 

2 2 
p i v i + 57 h i + 1L + ^w f + ¿w s = p 2 v 2 + h 2 * ’ÍL. 


2g 


(2.5) 


2 g. 


V Ap + Ah + 

g„ 


A\T 

— + . aw f + 4« s 


(2.6) 


Donde : 

V.- volumen específico medio del fluido (V - ) 

Multiplicando la Ec. (2.6) por p/ Al y considerando despreciables las 
perdidas de energía por trabajo externo, se tiene : 


AL H - + 


Av 


g c AL 


2g c AL 


+ p 


AW f 


(2.7) 


Considerando positiva la caída de presión en la dirección del flujo, 
se tiene: > 


AR - 

AL 


'4-#- + o- 


AV 


g AL 
c 


2g c AL 


+ P 


AW f 

AL 


( 2 . 8 ) 


A ésta ecuación se le acostumbra escribir en la forma siguiente: 


4EL) = ( AR-) + . ( AR s , , :_4P * 

Al T 1 AL e 1 AL ac + \ AL 


(2.9) 


Donde : 

( — ~ gradiente de presión total 


(-#-> 


¿L )' e ~ gradiente de Presión debido a la elevación. 


^~Al~ *ac = gradiente de presión debido a la aceleración. 
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( -t? ~ ) : - = gradiente de presión debido a la fricción. 

ZVL f 

2.2 Pérdidas de presión por fricción. 

Las pérdidas de presión por fricción, en conductos circulares de diámetro 
constante, han sido determinadas experimentalmente por varios investigado^ 
res. Los resultados de esta experimentación, utilizando tuberías de diver_ 
sos materiales, constituyen la base de las fórmulas que actualmente se - 
usan. 


2.2.1 Ecuación de Darcy. 

Darcy, Weisbach y otros, en 1857, dedujeron experimentalmente la siguien- 
te ecuación: 


Ap , = f P v 2 

( AL ‘f 2 g d 

c 


( 2 . 10 ) 


Una ecuación similar fue establecida posteriormente por Fanning, quien - 
obtuvo valores de f cuatro veces menores que los de Darcy. Esta dife- 
rencia se debe a el uso del radio hidráulico en lugar del diámetro de la 
tubería al formular su correlación. 


La ecuación de Darcy se usará generalmente, en estos apuntes, para calcu- 
lar las pérdidas de presión por fricción. 


2.2.2 Ecuación de Fanning. 

v 

La ecuación establecida por Fanning es : 

( .. Ap _) = f . h 

AL ' f 2g c 

Donde : 


R^ = radio hidráulico _ área de la sección transversal 
n perímetro mojado 


R^ = ( n d /4) / 7T d = d/4 


a 2 - ■) 
AL 1 f 


2 f v 2 p 


( 2 . 11 ) 
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2.2.3 Factor de fricción. 

El valor del factor de fricción (f > es función de la rugosidad de la tu- 
bería ( e ) y del numero de Reynolds (N ) , esto es : 


f ( •€ ', 


V 


El numero de Reynolds (adimensional) se define como : 

d v P 


N = 
Re 


P 


(2.13) 


Para calcular el valor de f , es necesario determinar el regimen de flujo. 
En conductos, los fluidos se mueven de acuerdo a cualquiera de los si- 
guientes regímenes de flujo : laminar o turbulento. El flujo laminar ocu- 
rre cuando las partículas de fluido se mueven en líneas rectas paralelas 
al eje del conducto. A velocidades mayores, las partículas se mueven de 
una manera caótica, formando vórtices y remolinos; en este caso el flujo 
es turbulento. Osborne Reynolds estableció experimentalmente un parámetro 
para determinar el regimen de flujo en tuberías.. A este parámetro (ecua- 
ción 2.13) , se le conoce como numero de Reynolds, El flujo laminar se pre 


senta cuando N < 2 300 
Re 


el flujo turbulento cuando N 


3100 


Para flujo laminar de una sola fase, el factor de fricción depende exclu 
sivamente del número de Reynolds y está dado por: 


•f - 


64 

Re 


(2.14) 


Para flujo turbulento (N Re — 3100) , el factor de fricción está dado por 
la ecuación de Colebrook y White : 




f = - 2 log ( 


+ 4^- i 


3.715 d 


(2.15) 


Se observa que para calcular f, en este caso, se requiere de un proceso 
iterativo. 

Basándose en la ecuación (2.15) Moody preparó un diagrama para determinar 
el factor de fricción en tuberías de rugosidad comercial (Fig. 2.2). 

En éste se nota lo siguiente : 

a) Para N Re <C 2300 (Flujo laminar) f depende exclusivamente del 
número de Reynolds . 


b) A partir de N =3100, se inicia la zona de transición. Dentro de 


ésta, f depende tanto de como de e /á (rugosidad relativa) 
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c) La zona francamente turbulenta se inicia a diferentes valores de N , 
dependiendo del valor de € /d. En esta zona f es independiente de Re 
N Re y varía únicamente con la rugosidad relativa. El valor de f puede 
obtenerse, para flujo turbulento, con: 

f = t - 2 log (í/ 3.715 d) )~ 2 (2.16) 


d) Cuando el flujo es crítico C2300 <N Re <3100) el factor de fricción 
se puede aproximar con la siguiente expresión: 

-O 3192 

f = 0.5675 N (2.17) 

Re 

Los valores de f , expresados en estas ultimas cuatro ecuaciones (2.14 a 
2.17) se utilizarán, junto con la ecuación de Darcy (2.10) , en el cálcu- 
lo de las perdidas de presión por fricción. 

Como se indicó, el valor de f, para flujo turbulento, es función también 
de la rugosidad e . Para tuberías comerciales e varía de 0.0006 a 0.0008 
pg. Para tuberías de producción comunmente se emplea un valor de e> 0.0006 
y para líneas superficiales de 0.0006 a 0.00075 pg. 


La siguiente ecuación permite obtener un valor de f bastante aproximado, 

cuando el regimen de flujo es turbulento (N >3100) 

Re 

(2.18) 
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Fig 2.2 Diagrama de Moody 



Coeficiente de fricción f 


8 




FLUJO DE LIQUIDOS POR TUBERIAS. 

Ecuación general en unidades prácticas. 

La ecuación (2.9), que expresa el gradiente de presión total, puie 
de escribirse en la forma siguiente, al considerar despreciable 
el efecto de la aceleración: 


¿P T = 4p e + 4p f 


La caída de presión por elevación, es: 
Ap ■* 0.443 y dh , 


Donde : Ap está en lb/pg , Ah en pies y V es la densidad 

relativa del líquido. r 

La pérdida de presión por fricción, en unidades prácticas, se ob- 
tiene con la ecuación de Darcy, de la manera siguiente: 

(3 3 , 

P f 2 g d' 

c 


p - 62 . 428 V r í 


/ EiÉS 

,.2 seg 


2 42 q f 2 

; , 

Sustituyendo (3.4) y (3.6) en (3.3) , se tiene: 

Api , l,S72768f (3 . 7 ) 

d 5 

Para emplear unidades prácticas se hacen las siguiente sustitu- 


d'5 (pies 5 ) « d 5 Cpg 5 ) 1 , pies 5 . 

,,5 1 5 ' 

12 pg 


3.1 


ftujo d<L Líquido ¿ pon. TubeAXcu 


rr' ( SÁ^Í y _ - r— V 5 ■■¿irffí f pies 3 s X , día 

q V seg } * W 5.6146 C ^ ■ - > 


•2 _ q 2 (5.6146) 2 


L' (pies) = L (millas) 5277 ) 

milla 


P f ) (^ ) 12 2 (-B2_) 

pie pg pir 


Sustituyendo de (3.8) a (3.11) en (3.7) se tiene: 

Ap - 1.572768 (5.6142) 2 (12) 5 (5277) f P r q 2 L 
f (12) 2 (86400) 2 d 5 


Ap f - 0.06056 f f r <Z 2 L 
, d 5 


Al sustituir (3.13) y (3.2) en ( 3.1) , se obtiene : 


^P T = 0.433 Y r ¿h + 0.06056 f * r 3 L 

d 5 


3. 2 Número de Reynolds 

Conviene también obtener una ecuación del número de Reynolds, en la 
que sus factores estén en unidades prácticas: 


Las sustituciones de unidades se hacen de la forma siguiente: 

... . 


d ' fJL ' 


3.2 


flujo dd Líquidos pon TuboJiXaÁ 


Sea: 


. 3 — „ .3 

/P.W.__ v _ _ fbl v „ pies 1 ,dia % 

q ( ség " q C día } 5 * 6142 C "^T^ ) 16400 W* 


(3.18) 



(3.19) 


Ib 

u' (-: m , - ■) = fi (cp) 0.00067197 ( _Í5 e 
” pie-seg 


(3.20) 


pie-seg cp 


d' (pies) = d(pg) 


pie 


) 


12 pg 

Sustituyendo de (3.18) a (3.21) en (3.17) se tiene: 


(3.21) 


Re 


= 92.2 qjr 


d ¡Ji 


(3.22) 


3.3 Eficiencia de flujo. 1 

Es obvio que la rugosidad de las tuberías dependerá del proceso em- 
pleado en su fabricación, su grado y tamaño. Aun las tuberías nuevas 
y con mayor razón las almacenadas, mostrarán valores aleatorios. en 
su rugosidad. Los efectos de la corrosión, erosión e incrustamientos 
que ocurren en las tuberías en operación, también afectan las condi- 
ciones de flujo. Por lo anterior los gastos calculados mediante las 
correlaciones raramente concuerdan con los medidos. 


Para compensar esta imprecisión, generalmente se introduce en los - 
cálculos un factor de eficiencia E. Este factor se define como la 
fracción (o porciento) del gasto total calculado al manejado real- 
ícente en una tubería. Cuando se carece de datos de campo, la selec- 
ción de E se basa en la experiencia,* sin embargo, un valor de 0.90 
es apropiado para los cálculos en el diseño de tuberías. 

Para calcular el gasto real de una línea, su gasto teórico se multi- 
plica por el factor E. Para corregir las perdidas de presión calcu- 
ladas , estas se dividen por la raiz cuadrada de E. Este procedimien- 
to, tiene el mismo efecto que cambiar el factor de fricción por un 
nuevo valor f/E^. 

Por lo tanto la Ec. (3.14), de acuerdo, con lo expuesto , queda : 


Ap_ = 0.433 

T 


Ah + 0.06056 f ^r q L 


_2 . 5 
E d 


(3.23) 
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3.4 Aplicación de la ecuación de flujo en el análisis y diseño de tube- 
rías que conducen líquidos., 

La^ecuáción (3.23) permite calóular con aproximación aceptóle la - 
caída de presión en tuberías que conducen líquidos; ó bien, dada la 
caída de presión disponible, calcular el gasto o el diámetro. La 
ecuación es aplicable tanto al casó: de tuberías superficiales que 
conducen aceite (oleoductos) , como al de pozos inyectores de agua. 

3.4.1.- Ejemplo de cálculo de 


Calcular la caída de presión por fricción en una tubería de 3000 
pies de largo, 3.937 pg de diámetro interno y con una rugosidad ( e ) de 
0.0006 pg, donde fluye aceite de densidad relativa de 0.9 y /u =46 cp, 
si el gasto es de: a) 2560 bl/dia; b) 12 800 bl/dia. 

Solución : 


a) 


(3000 pies) 1 milla 

5277 pies 


- 0.5685 millas 


Re 


92-2 g y _r (92.2) (2560) (0.9) _ . 

<3 /< (3.937) (46) 1173 


Ya que N <12300, el flujo es laminar y f se calcula con la ecuación 
(2.14) : Ke 


f = 


64 

N 

Re 


64 

1173 


- 0.05456 


Ap f ^ 0.06056 fjr gf L 
d 5 


An (0.06056) (0.05456) (0.9) (2560)^ (0.5685) 

kf 5 

(3.937) 3 


¿P f = 11.7 lb/pg 


b) Rara el gasto de 12 800 b 1/día se tiene: 
(92.2) (12 800) (0.9) 




Re 


(3.937) (46) 


- 5865 


ya que N >.3100, el flujo es turbulento. El factor de fricción, según 
la ecuación (2.15), se puede obtener con el siguiente procedimiento. itera- 
tivo : 
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f = 


-2 log 


3.715 d 


2.514 


s Re 


-2 


(3.24) 


En donde f s es el valor supuesto y ^es el calculado. Esta ecuación se ~ 
aplica a partir de un valor supuesto de f (f ) obteniéndose un valor más 
aproximado f c . Si j f c - f g j ^tolerancia, S f c es el factor de fricción 

buscado; si no, se hace f g = f c y se repite el procedimiento. 

Para este ejemplo: 


f c = - 2 


[0.0006 

, 2514 

|1 

l (3.715) (3.937) 

5865 '¿T 

s 

/] 


-2 log |ó 


.000041 + 0.0004286 


V f S 


1-2 


Iniciando con f =0.03, se obtienen los resultados que se muestran a con 
t limación : 

Iteración f 


0.03700 

0.03576 

0.03596 

0.03592 

0.03593 


Por lo que: 


Ap = 0.06056 f *r q L 

,5 


A-p - 


(0,06056) (0.03593) (0.9) (12 800) (0.5685) 
(3 . 937) 5 


Ap f = 192.8 lb/pg 


3.4.2 Diseño de tuberías. 


La ecuación (3.23) puede aplicarse también para obtener el diámetro para 
un gasto y caída de presión dados. El procedimiento de solución es itera- 
tivo, ya que el numero de Reynolds (y por lo tanto el factor de fricción ) 
es función dél diámetro. V* 
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De la ecuación (3.23), despejando d se obtiene: 



Como se indico el procedimiento iterativo que se usa para resolver estas 
ecuaciones es el de sustituciones sucesivas; o sea, para una iteración 
dada se supone como aproximación de la incógnita el valor calculado en la 
iteración anterior. Las iteraciones se continúan, hasta que el cambio 
de la incógnita en una iteración, sea menor que la tolerancia estableci- 


Se observa que es necesario realizar iteraciones al calcular el diámetro 
o el gasto y al calcular el factor de fricción cuando el flujo es turbu- 
lento. De esto se desprende la conveniencia de usar calculadoras programa- 
bles en la solución de problemas de diseño y análisis de tuberías horizon 
tales, verticales o inclinadas. 


3.4.3 Predicción del comportamiento de pozos inyectores de agua. 

La Fig. 3.1 es un diagrama simplificado de un sistema de inyección de - 
agua. Los elementos del sistema son: 

a) La estación de bombeo 

b) La línea superficial de distribución 

c) La tubería vertical de inyección 

d) El yacimiento 

En la misma figura se distinguen los nodos o puntos de unión de los ele- 
mentos, en los que se tienen las presiones siguientes : 

a) Pp*~ Presión en la estación de bombeo. 

b) p-, .- Presión en la boca del pozo 

tn 

P wf.- Presión de fondo fluyendo (inyectando) 

d) p .- Presión estática (de fondo cerrado) . 
r ws 
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Estación 
de bombeo 


Pp 
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La predicción del comportamiento de un pozo inyector consiste esencialmen 
te en la^determinacion de su capacidad de admisión (gasto de inyec ción ) ~ 
en función de la presión disponible en la superficie y las características 
de las tuberías y del yacimiento. Los datos requeridos para resolver el 
problema planteado son: 

a) Presión de inyección en la estación de bombeo. 

b) Longitud, diámetro y cambio de nivel de la tubería superficial de 
distribución. 

c) Longitud y diámetro de la tubería vertical de inyección. 

d) Presión estática e índice de inyectividad (II) del pozo. 

e) Características del fluido inyectado. 


El procedimiento de cálculo consiste en : 

a) Suponer un gasto de inyección . 

b) Con la presión de ^bombeo y las características de la línea superficial, 
calcular la presión en la boca del pozo para el gasto supuesto. Para 
este cálculo, se usa la Ec. (3.23) . 

c) Con la presión en la^boca calculada en el paso anterior y el mismo gas 
to, obtener la presión de fondo inyectado, (p ) . para este cálculo 

se usa la misma Ec. (3.23) , teniendo en cuenta W que el desniveldh (pies) 
corresponde al valor de la longitud de la tubería vertical de inyec- 
ción L, que debe estar en millas. 

d) Con el índice de inyectividad del pozo, el gasto supuesto, y la presión 
de fondo calculada en el paso anterior, obtener la presión estática 


q 

II 


(3.27) 


e) Comparar el valor de la presión estática calculada, con el valor real 
de dicha presión. Si estos valores coinciden, dentro de la tolerancia 
fijada, el gasto supuesto en el paso (a es el correcto . En caso con- 
trario, suponer otro gasto y repetir el procedimiento, hasta obtener 
la aproximación deseada. 

Para simplificar los cálculos es conveniente graficar los valores de 
los gastos supuestos contra los de la presión estática calculada. Co- 
mo se muestra en la Fig. 3.2. 
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FIG. 3*2.- DETERMINACION GRAFICA DEL GASTO CORRESPONDIENTE A UNA 
PRESION ESTATICA DADA. 


El procedimiento descrito permite obtener la presión en cada nodo, en - 
función de las condiciones de flujo que prevalecen en los elementos del 
sistema. Por esta razón el procedimiento se denomina como análisis nodal. 

Con relación al comportamiento de un pozo inyector de aguaces importan- 
te considerar que su presión estática generalmente aumentara con el tiem- 
po, debido a la propia inyección. Por otra parte la formación puede sufrir 
daño, por las impurezas del agua o su incompatibilidad con la formación. 
Además la movilidad del agua en la zona invadida, puede ser menor que la 
movilidad de los fluidos desplazados delante del frente de invasión. El 
efecto de estas situaciones es la disminución del índice de inyectividad. 
Por lo tanto es conveniente determinar la capacidad de inyección del po- 
zo considerando como variable la presión estática y el índice de inyecti- 
vidad. 

Es evidente que un análisis detallado permitirá determinar el efecto que 
sobre el gasto, tendrá el cambio de las tuberías y la presión de inyec- 
ción. 


3.4.4 Ejemplo de aplicación del análisis nodal. 

Se desea determinar: a).- El gasto de inyección de un pozo; y b) El efec- 
to del cambio de la presión de bombeo y de las tuberías, sobre su capaci- 
dad de admisión. Los datos del sistema son: 


3. 9 


Flujo de Líquidos pon. Tube/Uaó 

Para la línea superficial: 

= 1.00 L = 2 millas ju= 1 cp 

E = 0.90 d - 4 pg Ah = + 20 pies € = 0.0006 pg 

Para la tubería vertical de inyección: 

y = 1.00 L = 2.273 millas /“ = 0.8 cp 

rw 

E =0.90 d = 2.441 pg Ah = - 12000 pies 

e= 0.0006 pg 

2 

Otros datos son: p^ = 3100 lb/pg 

2 2 

II = 3.521 bl/día/lb/pg P ws “ 4402 lb//pg 

Solución : 

a) Siguiendo el procedimiento descrito se obtiene los valores de p wg 
mostrados en la tabla siguiente: 

TABLA 3.1.- CALCULO DEL GASTO DE INYECCION 


GASTO SUPUESTO 




(bl/día) 

P th 

P wf 2 

P ws 



(lb/pg ) 


6000 

2997.8 

7011.3 

5307.2 

8000 

2932.5 

6084.4 

3812.3 

10000 

2851.0 

4913.1 

2073.0 

7250 

2959.0 

6460.7 

4401.6 

En la Fig. 3.3 

se muestran los 

valores anteriores 

graficados. Se 

va que el gasto 

de inyección es 

de 7250 bl/día. 



3.10 


Flujo d<¿ Líquido^ pon. TuboAMU 


b) En la Tabla 3.2 se presentan los valores de los gastos calculados, si- 
guiendo el procedimiento enunciado, al suponer tuberías, de diferentes 
diámetros y presiones de bombeo variables. Los resultados se muestran 
en las Figs. (3.4), (3.5) y (3.6). 


TABLA 3.2. EFECTO DE LAS TUBERIAS Y LA PRESION DE BOMBEO SOBRE EL 



GASTO DE INYECCION. 




L.D. 

q 

T.P. 

q 

P P 2 

q 

(pg) 

(bl/día) 

(pg) 

(bl/día) 

(lb/pg > 

(bl/día ) 


2 

4690 

2 

5105 

2000 

5710 

3 

6765 

3 

9422 

3100 

7250 

4 

7250 

4 

11400 

4000 

8375 

5 

7365 

5 

12050 

5000 

9530 

6 

7405 





Es obvio que la selección 

de las 

tuberías debe 

basarse en 

un análisis eco- 

nómico , 

, en el que se comparen los incrementos 

en el gastO; 

, al instalar tu- 

herías 

de mayor diámetro, 

con la 

inversión adicional que es necesario rea- 


lizar. 


3.4.5 Perdidas por fricción en las conexiones. 

Los codos, las válvulas y las conexiones, así como otros tipos de adita- 
mentos, incrementan las pérdidas de presión por fricción en el sistemé. 

En el mejor de los casos, éstas pérdidas adicionales de presión solo pue- 
den ser estimadas en forma aproximada. Ellas deben de incluirse en el ana 
íisis de un sistema, considerando que cada conexión es sustituida por una 
longitud equivalente de tubería recta (Le) , la cual producirá la misma 
pérdida de presión por fricción que la conexión real. La longitud equiva- 
lente de cada conexión se agrega a la longitud axial de la tubería (L) , 
antes de calcular la pérdida de presión total en el sistema. Los valores 
aproximados de Le, para las válvulas comunes y otras conexiones, para flu- 
jo turbulento, pueden obtenerse mediante el nomograma mostrado en la Fig. 
(3.7). Por ejemplo un codo de 2 pg. , en ángulo rectores equivalente a 12 

pies de una tubería recta. 
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Fig 3.4 Relación entre el gasto de inyección y el diámetro de la 
línea de descarga 



Fig 3.5 Variación del gasto de inyección ai considerar tuberías 
de produccio'n de diámetros diferentes. 



Fig 3.6 Relación entre el gasto de inyección y la presión de 
bombeo 
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4. FLUJO DE GAS POR TUBERIAS Y ESTRANGULADORES . 

4. I. - Ecuación general en unidades prácticas. 

La perdida de presión por fricción, según la ecuación de Darcyj^ es: 

_ 1 ÍL ¿ k ‘ <4 i) 

áP f ~ 2 a d' ( ’ 

c 

Para expresar esta ecuación en unidades practicas, se haran las substi- 
tuciones siguientes: 


De la EC. (1. 36) : 


2.7044 p V a 
Z (T + 460) 


v = q'/A = 


Bg, según se vio (Ec. (1.34)), es: 


0.02825 Z (T'+ 460) 


Sustituyendo (4.4) en (4.3) y elevando al cuadrado: 


2 q' 2 Z 2 (T + 460) 2 

v — 4~l2 

772.93 d p 


Sustituyendo (4.2) y (4.5) en (4.1) y simplificando: 

f V* g' 2 Z L* (T + 460) 

APf = 5 - 

18389 d' p 


Hasta aquí los términos de la ecuación, excepto p, están en unidades ini- 
ciales, o sea: 

p' (lb/pie 2 ) , q* (pies 3 /seg) , ¿(pies) , d*(pies) * 
p ( lb/pg 2 ) , T ( °F) 

Las sustituciones, por unidades prácticas, se hacen en la forma siguien 
te : 

144(22L> (4.7,1 

pie P9 P ie 
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d ,s (pies = d Cpg 5 ) < pg5^"“ * 

12 pg 


(4.8) 


naf 2 - 1 


1 - I-S^l 


2 v 2 
86400 seg 


(4.9) 


Además: p 
Donde : 


L' (pies) * L (millas) 5277 ( > 

P 1 + P 2 


milla 
y 4p = p x - p 2 


( 4 . 10 ) 


ss presión a la entrada del ducto 
p^ =s presión a la salida del ducto 
Haciendo las substituciones anteriores en la Ec. (4.6) se obtiene: 


f v q Z L ( T + 460) 
. 2 2 . Y q 

(P 1 “ P 2 ) f 


2744 2 d 5 


(4.11) 


Esta ecuación permite obtener las perdidas de presión por fricción. 
Ahora se requiere una expresión para obtener las perdidas de presión 
por elevación ( A P )• Para ello se procede de la siguiente manera: 


AP a ( 


Jf-> = f ( 


Ib 


. 2 


pg” Pie 

Donde h es el cambio de altura. 


144 % nrT 2 
pg 


(-Eí|_) h (pies) (4.12) 


Como: 


= e 


2.7044 p \r q 
Z (T + 460) 


(1.36) 


y : 

AP e - Pi - *>2 


Sustituyendo (1.36) y (4.13) en (4.12) : 


(p l ‘ Ve 


0.01878 P V CT h 
Z (T + 460) 


(4.13) 


(4.14) 


Multiplicando ambos miembros de la igualdad por p (p^ + ¥ 

despejando p 2 _ p 2 
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2 2 
( * i - p 2 > e 


-2 

0.Q3756 p h 

Z (T + 460) 


<4.15) 


Sumando (4.11) y (4.15) se tiene: 


Pr - Po = 


f V q q 2 Z L (T + 460) + 0.03756 p 2 V q h 

2744 2 d 5 Z <T + 460) 


Con ésta ecuación pueden calcularse los valores de p^ o p^. Sus. unida- 
des son: 

2 3 

p (lb/pg ), q (pies^ a c.s/día) , L (millas), d (pg) 
h (pies) , T ( °F) 


Despejando q de la ecuación anterior se obtiene: 


*= 2744 ( ^ ) 0-5 d 2.5 


P 1 * P 2 - (0 - Q375 V q h P 2 ^ <£ + 460 >) 


7 . LZ (T + 460) 


0.5 


Asimismo se tiene: 
2 


d = 


q f yv/ L Z ( T + 460) 


2744 2 ( p* - p 2 2 - (0.0375 y h P 2 /Z (T + 460))) 


4.2 Numero de Reynolds . 

d 1 v* p' 


(4.19) 


Según se ha visto: 

_ 2.7044 p ? a ._ 

r Z (T + 460) 


4 q' B ^ (4) (0.02825) q 1 Z (T + 460) 


7t d 


* d ' 2 p 


(4.20) 

(4.21) 


Sustituyendo (4. 20) y (4 .21) en (4. 19) 

V <3 


N„ * 0.09727 q< F( 3 
Re d' ju’ 


(4.22) 


(4.16) 


(4.17) 


2 (4.18) 
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Procediendo ahora al cambio de unidades: 


q « (JBÍSSÍ ) . q ( ) 

4 1 c 0 n ' 4 v dia 


seg 


— -i— (-JÜ* .) 

86400 s seg 


(4.23) 


d* (pies) * d (pg) ~ (- £ ~~ ) 


(4.24) 


ib 


Ib 


M'(- 


pie - seg 


) = }i (cp) 0.00067197 ( 


pie- seg cp 


(4.25) 


Sustituyendo estas tres últimas igualdades en (4.22) y simplificando: 

N - 0.020105 § ■■ -S - - (4.26) 

Re d « 

Mg 

4.3.- OBSERVACIONES: 

4.3.1. Como la viscosidad del gas es baja y generalmente se manejan gas_ 
tos altos, el regimen de flujo es normalmente turbulento, por lo que 
el valor de f, según la ecuación (2.16) , es: 

f = ( 2 log (3.715 d/e) ) “ 2 (4.27) 


4.3.2. Al emplear las ecuaciones (4.16), (4.17) y (4.18) se debe tener 

cuidado en la determinación del signo de h. Para esto basta recordar 
que la diferencia de alturas entre la descarga y la entrada se debe re 
ferir a un nivel base y se considera el sentido positivo hacia arriba. 
Así para un pozo inyector, en el que la entrada esta en la superficie 
y la descarga en el fondo, a una profundidad D, la altura de la entra- 
da respecto a la superficie es cero y la de la descarga referida al mijs 

mo nivel es -D, de tal manera h * h D - h e = - D - 0 - - D. Lo mismo 

se obtiene si el nivel de referencia es la profundidad D. En este caso 

h D = 0yh/=D, o sea: h = h D - hg = 0 - D = - D. 

4.3.3. Para resolver la ecuación (4.16) se requiere del uso de un pro^ 
cedimiento iterativo, ya que p, z y f dependen de la presión y f ademas 
depende del gasto. El procedimiento iterativo recomendado es el de apro 
ximaciones sucesivas, mencionado en el inciso 3.4.2. 

4.3.4. Como la distribución de presiones en un gasoducto no es lineal, 
se recomienda el uso de la siguiente ecuación para obtener la presión 
media: 

2 

P - — 


P 1 + P 2 


(4.28) 
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4.4.- Predicción del comportamiento de un pozo productor de gas. 

Sr^» 9rama simplificado del sistema de producción de un pozo produc- 

v LfcIraJT “ " U - Stra “ U Fig - 4 ' X - elementos del síteme 
y las caídas de presión que en ellos se presentan son: 


ELEMENTOS 
1*- Yacimiento 

2 . - Tubería de producción 

3 . - Estrangulador 

4. ~ Linea de descarga 


PERDIDAS DE PRESION 

D — n 

WS *wf 

p vf - P th 




FIO- 4.1.- SISTEIíA SIMPLIFICADO DE FLUJO DE UN POZO PRODUCTOR DE GAS 
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La capacidad de producción del sistema depende de la presión estática 
del pozo, las características del yacimiento, las tuberías, el estran- 
gulador, la presión de separación y las propiedades del gas producido. 

Con la información anterior se puede obtener el ritmo de producción me_ 
diante el procedimiento siguiente: 

a) Suponga un gasto 

b) Con la presión, estática y la ecuación que representa el flujo en el 
yacimiento,' obtenga, para el gasto supuesto, el valor de la presión 
de fondo fluyendo. 

En estos apuntes se considerará que el flujo en el yacimiento puede 
quedar representado por la ecuación siguiente: 

q .= c (P^s - 4>" ¿ ' (4-29) 

Despejando se tiene: 

p , =. fr? - ( a ) 1/n (4.30) 

wf y ws c 

c) Con la presión de fondo fluyendo y el mismo gasto, obtenga la pre- 
sión en la boca del pozo, p^» Para. este cálculo se usa la Ec. 4.16. 

El valor obtenido representa la presión antes del estrangulador . 

d) La secuencia de cálculo se inicia ahora a partir de la presión de 
separación , para obtener la presión en la boca del ,; pozo corriente 
abajo del estrangulador (p ), necesaria para transportar el gasto 
supuesto a través de la línea de descarga. 

e) Repita el procedimiento suponiendo diferentes gastos y grafique los 
Valores de p^ y de p^ contra el gasto, como se indica, en la Fig. 4.2 


En dicha figura se observa que cuando p = p^ se tiene el gasto máxi- 
mo, correspondiente al flujo sin estrangulador. Para obtener gastos 
inferiores al anterior se requiere un estrangulador en el cabezal del 
pozo. El tamaño del estrangulador se puede calcular con las ecuaciones 
que se presentan en el inciso 4.5. 

f) En la misma forma, calcule los gastos máximos que pueden obtenerse - 
para diferentes presiones estáticas. De esta manera se determina la 
relación entre estas variables, lo que constituye en sí la predicción 
del comportamiento del pozo. 
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FIG. 4.2.- RELACIONES CALCULADAS DE LAS PRESIONES A BOCA DE POZO PARA 
DIFERENTES GASTOS . 


4.4.1. Ejemplo de predicción del comportamiento de ion pozo productor 
de gas seco. 

Aplique el procedimiento descrito para las condiciones de flujo s iguien 
tes : — 

Presiones estáticas de 2000, 1500, 1000 y 500 lb/pg 2 

. ^ ' 2 
Presión de separación de 50 lb/pg 

Diámetro de la tubería de producción : 2 . 441 pg 

Diámetro de la línea de descarga: 2.991 pg 

Longitud de la tubería de producción: 6320 pies = 1.1976 millas 
Longitud de la línea de descarga: 1.865 millas. 

Desnivel del separador: 300 pies (flujo ascendente) 

Densidad relativa del gas: 0.70 (aire ' ss. 1.00) 

Temperatura media del fluido en la tubería de producción: 140°F 
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Temperatura media del fluido en la línea de descarga: 60 °F 

El comportamiento del flujo en el yacimiento está dado por la Ec. (4.30) 
donde: C = 16280 y n - 0.681 

La rugosidad de las tuberías es de 0.0007 y se considera E = 1.00 

Al seguir el procedimiento expuesto se obtuvieron los siguientes valo 
res de P th / partiendo de las presiones estáticas consideradas: ”” 


p = 2000 lb/pg 2 

ws ■ 


p = 1500 
ws 

q (pies 3 /dia) 

P wf 

p th 


q 

P th 

18 x 10 6 

1992.6 

472 


12 x ib 6 

650 

16 x 10 6 

1993.8 

886 


10 x 10 é 

979 

12 x 10 6 

1995.9 

1301 


8 x I0 6 

1052 


P ws - 1000 lb/pg 2 

g 

P th 

2 x !0 6 

, 834 

6 x 10 6 

642 

8 x 10 6 

401 


** 


p = 500 lb/pg 2 

WS 

q 

p th. 

1 x ío 6 

420 

2 x io 6 

385 

3 x lo 6 

317 

4 x 10 6 

185 


Para la línea de descarga se obtuvo el comportamiento siguiente : 


r — 

q 

- 

2 x io 6 

153 

4 x 10 6 

292 

10 x lo 6 

698 

12 x 10 6 

828 

16 x 10 6 

1081 


Al graficar los valores tabulados anteriores se obtuvo la Fig. 4.3. 

De esta, con los valores correspondientes al gasto máximo sé obtuvo la 
Fig. 4.4, que representa el comportamiento del pozo. 
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q g x 10"" s (p¡esVdía) 


Fig 4.3 Gastos obtenibles, a diferentes presiones estáticas, en fun 
ción de las presiones en el cabezal del pozo. 



Fig 4.4 Relación entre el gasto máximo y la presión estática 
del pozo. 







Fig 4.5 Distribución de las presiones en los nodos , en función 
del gasto, para una p ajS s 500 lb/pg 2 



Fig 4.6 Relación de calores específicos en función de ía tem- 
peratura y la densidad relativa 
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La aplicación del procedimiento descrito permite analizar el efecto so^ 
bre la producción al variar las tuberías de producción y de descarga, 
así como el de la presión de separación. 

La distribución de presiones en los nodos, para una presión estática 
dada, en función del gasto, se muestra en la Fig. 4.5. 


4.5.- Flujo de gas a través de éstranguladores 


( 2 ) 


Los estranguladores que se instalan en la boca del pozo, para controlar la 
producción, están basados en el principio de flujo crítico. 

El flujo crítico es un fenómeno de flujo definido por el flujo de gases 
comprensibles , en la sección de estrangulamiento de una restricción, 
cuando su velocidad es sónica (velocidad del sonido en el fluido) o el 
numero Mach es uno. El numero Mach es una relación adimensional dada 
por el cociente de la velocidad del fluido entre la velocidad del soni- 
do . 

Existe flujo subcrítico si: 


P 2 / P x > (2/K + 1). 


K/(K - 1) 


(4.31) 


Se tiene flujo crítico cuando: 


(4.32) 


Donde : 


calor específico a presión constante 
~ - calor específico a volumen constante 


El valor de K puede obtenerse de la Fig. 4.6 

Las ecuaciones siguientes permite calcular el diámetro del estrangula— 
dor. 


(4.33 


C A - 



■a • 



■ 

155500p 1 

64.34 K 

P, 2/K 

) ' 

p (k+D/k 

r A . \ 

0.5 

y (T + 460) (K-l) 

r.g 

( £>i J 



Donde: 

2 

A * área del estran guiador, pg 
q = coeficiente de descarga del estrangulador. 
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T . = temperatura en el estrangulador °F 

p^ ~ presión corriente arriba del estrangulador (p^) 

P 2 = presión corriente abajo del estrangulador (p^) 

El diámetro del estrangulador puede obtenerse con la siguiente ecuación, 
ajustada de la correlación establecida por Cook^ 2 ^. 

d 1.75105 + 932.334 CA - 2 9372. 7(CA) 2 

C 3 4 

+ 397972 (CA)' - 1510615 (CA) (4.34) 


Donde: 

d = diámetro del estrangulador en 64 avos de pg. 
c 

Si existe flujo crítico a través del estrangulador , el diámetro se pue_ 
de calcular con las mismas ecuaciones (4. 33). y (4.34); pero en lugar del 
valor de p^/p^ sé usará el de (2/K + 1) K//K 

Cuando se tiene flujo crítico a través del estrangulador/ la presión co^ 
mente arriba es independiente de la presión que prevalece en el siste- 
ma de recolección (línea de descarga, separadores, etc.) Se infiere que 
el control de la producción se logrará cuando las variaciones de pre- 
sión en el sistema de recolección no se reflejen en la formación pro- 
ductora, provocando fluctuaciones en la producción. Esta situación pre_ 
valecerá al usar un estrangulador que permita obtener la producción de 
seada bajo condiciones de flujo crítico. 

Ejemplo: Calcule el diámetro del estrangulador si: 

q « 80 000 pies 3 a c.s. T = 150°F 

9 

y — 0.6 K = 1.31 

* 9 

7 2 

p x = 700 lb/pg P 2 « 600 lb/pg 

Solución : 

( —1 . )K/(K-1) = 0#545 

K + 1 

P 2 / = 0.857 

Como p 2 / P L > (2 / K + 1 ) K/(K ” 1) , se usará el valor de p 2 / P x en 
la Ec. (4.33 ) , ya que el flujo es subcrítico. 
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CA 


80 000 


(155500) 700 


(64.34) (1. 31) 

(0.6) (610) (0.31) 


(0.857) 


1.527 


(0.857) 




0.5 


2 

CA * 0.005068 pg 

Sustituyendo este valor en la Ec. (4.34 ) , se obtiene: 

d = 5.77/64 pg -■ 6/64 pg. 
c 
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5.- FLUJO MULTIFASICO EN TUBERIAS, 


5.1. - Conceptos y ecuaciones fundamentales. 

5.1.1. Patrones de flujo.- Es evidente que al fluir dos fases simultánea 
mente, lo pueden hacer en formas diversas. Cada una de estas formas con£ 
tituye un patrón de flujo. La distribución relativa de una fase con res- 
pecto a la otra, se conoce como patrón o tipo de flujo. 

En el flujo vertical se han identificado patrones de flujo burbuja, tipo 
bache y tipo niebla, así como el correspondiente a la transición entre 
estos dos últimos tipos de flujo. 

En la Fig. 5.1 se presenta una clasificación general de los patrones de 
flujo observados durante el trabajo experimental realizado por Beggs^* 

5.1.2. Colgamiento (H L ) .- Se define como la relación entre el volumen de 
líquido existente en una sección de tubería a las condiciones de flujo, 
entre el volumen de la sección aludida. Esta relación de volúmenes depen 
de de la cantidad de líquido y gas que fluyen simultáneamente en la tu- 
bería. Generalmente la velocidad con que fluye el gas es diferente de la 
velocidad con que fluye el líquido, existiendo un "resbalamiento" de una 
de las fases. 

El término "resbalamiento" se usa para describir el fenómeno natural del 
flujo a mayor velocidad de una de las dos fases. Las causas del .resbala- 
miento son diversas. La resistencia al flujo por fricción es mucho menor 
en la fase gaseosa que en la fase líquida. La diferencia en compresibili 
dades entre el gas y el líquido hace que el gas en expansión viaje a ma- 
yor velocidad que el líquido. El resbalamiento también es promovido por 
la segregación gravitacional . Las fuerzas gravitacionales originan que 
el líquido se mueva a mayor velocidad que el gas, cuando el flujo des- 
cendente; pero ocurre lo contrario en el flujo ascendente. 

Para calcular las pérdidas de presión por elevación (carga hidrostática) , 
es necesario predecir con precisión el colgamiento (H^) considerando el 
resbalamiento entre las fases. 


Existen varias correlaciones para obtener el colgamiento del líquido. La 
correlación mas general es la de Mukherjee y Brill,^ ' obtenida a partir 
de mas de 1500 mediciones para flujo con ángulos de inclinación de 0 a 
+ 90°. Los resultados de este trabajo son los que se usarán en estos - 
apuntes z 


La ecuación, 

H = exp 
L 


establecida por Mukherjee y Brill, es: 
2 2 

(C, + C sen 0 + C_ sen 0 + C . N ) 

1 2 3 4 L 


K CS 

_£L 


C6 

tC 

Lv 


(5.1) 


5.1 


Flujo 

segregado 


Flujo 

intermitente 


Viujo 

distribuido 




Estratificado 


Ondulado 


Anular 


Tapón 


Bache 


Burbuja 


Niebla 


Fig 5.1 Patrones de flujo observados por Beggs en flujo 
horizontal 

Se tendrá flujo descendente estratificado cuando: 

. 1() (0.321 - 0.017 - 4.267 sen e - 2.972. N L 

LV — 

2 2 

- 0.033 (loq Ni - 3.925 sen ©) 

g 7 
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N_ *■ numero de la viscosidad del líquido * 0,15726 ¡jl 11/ p o- 3 ) °’ 25 
" L L 


*■ numero de la velocidad del líquido » 1.938 v C p ' / <r ) 


N gv * numero de la velocidad del gas = 1.938 v sg ( p a ) 


Las unidades usadas en los términos de estos números adimensionales son: 
a (c.p) , P (lb/pie 3 ) , <r (dinas/cm) v( pies/seg) 


Los coeficientes de la Ec. (5.1) para flujo ascendente y horizontal, de£ 
cendente, y para flujo estratificado descendente, se muestran en la ta- 
bla siguiente: 


Dirección Tipo de 
del Flujo Flujo 




n 

c: 

2 

3 

4 

0.12988 

-0.11979 

2.34323 

4.80814 

4.17158 

56.26227 

0.78981 

0.55163 

15.51921 



Otro concepto que se usa con frecuencia en los cálculos de gradientes pa 
ra flujo multifásico, es el colgamiento sin resbalamiento ( A ) . Se defi- 
ne en la misma forma que H L ; pero se calcula a partir de las condiciones 
de flujo (p y T) existentes , considerando las producciones obtenidas en 
la superficie (q y R) , esto es : 


q L + q g 


qn (R - V W 
5.615 (q B + q B V 
.o o *w w ■ 


Donde q = gasto a condiciones de escurrimiento . 


5.1.3.- Velocidades superficiales.- Es la velocidad que tendría cualquie 
ra de las fases si ocupara toda la tubería. Se define por las expresio- 
nes siguientes: 


v gL (pies/seg) 


v (pies/seg) 
. sg 


0.002122 o (R - R ) B 
. • • “o s g 
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Donde; 


A .- Area de la secciSn transversal de la tubería. 
P 


q T + q n 

h 2 _ 


m 


■ v _ + v 

sL sg 


(5.5) 


De estas ecuaciones se observa que; 

V sL 


.1 : 


(5.6) 


m 


5.1.4.- Velocidades reales.- Aplicando el concepto de colgamiento, se 
pueden obtener las velocidades reales de cada fase : 


q ¿ 


q L 


A H 
P L 


(5.7) 






A (1-H ) (1-H ) 

p L L 


(5.8) 


5.1.5.- Densidad de la mezcla de fluidos 


La densidad real de la mezcla de fluidos se obtiene a partir del colgá- 
is + P (1 - H ) (5.9) 

1 q L 


miento (H ) con : 

If 


Pm ” f*L "L lg **L' 


Algunos autores calculan la densidad de la mezcla sin considerar el res^ 
balamiento entre las fases/ esto es: 


P = P T A + P (1 - A ) 

Vns V L v g 


(5.10) 


También puede obtenerse esta densidad a partir de la expresión siguien 
te ; 


e _ £L_ 

ns V 


(5.11) 


Donde ; 

M - masa de la mezcla a c. esc. por barril de aceite producido a c.s. 
(Dbm/blo a c.s.) 
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V m “ volumen de la mezcla a c. esc, por barril de aceite producido 
pies^ a c. esc., /b lo a c.s. 


loas valores de M y V se obtienen con las ecuaciones siguientes: 


M = M + M + M 
o g w 


(5 


. 3 


«O • r„ > . >•<*» ¿SÍ- . 


(lbw/pie J ) 


P ie w 


M = 350.5 y 


(5 


Ib /pie ) 

M y -L_a a_ 

g g (lb a /pie^) 


Ib 


x 0.0764 ó 


, 3 
pie 

a 


pies a c.s 

) X R (~¡ r3E ) 

bl a c.s 


M 0.0764 y R 
g * g 


(Ib / pie ) 

y Sí i-J i— 

w , / . 3 , 

( Ib / pie ) 


(5 


. 3 

Ib pie 

x 62.428 ( ) x 5.615 ( -- --- 

J bl 

pie w 


• bl 

x WOR ) 


M = 350.5 y WOR 
w w 


(5 


. M = 350.5 V + 0.0764 R y + 350.5 


WOR 


(5 


Cálculo de V (- 
m 


pies^ + g + w a c. esc. ) 


bl a c.s. 


V * 5.615 B 

mo o 


v = (R - R > ( pies 3 gl a e. esc, a c.s ) B c Pies 3 c. esc. 


mg 


blo a c.s. 


g . 3 

pie a c.s. 


a c.s. 

. 12 ) 

.13) 

.14) 

.15) 

.16) 
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. 3 

V - WOR (Z±2Í — i 5.615 c “TT— 
mw . , b/w 

Di 


) b* (• 


pies 


w a c. esc.. 


pie a c. 


V “ 5.61 B + (R - R ) B + 5.615 B WOR 
m o s g w 


Sustituyendo (5.161 y (5.17) en (5.11), se obtiene 

350.5 ( y + v WOR) +0.0764 R y 
r ro Y w , g 

^ns “ 5.615 CB q + B w WOR ) + (R-Rg) B g 

5.1.6.- Gasto de masa.- Se define por la siguiente expresión: 
Ib de líquido y gas 


m segundo 

Puede obtenerse con cualquiera de las ecuaciones siguientes: 
q M 


m 86400 


w = w + w + w 
m o w g 


w = p IB/ 15391 
o o o o 


W = P 


W %í w 


B / 15391 


o q (R - R ) B / 86400 
■ . o s g 


5.1.7.- Viscosidad de la mezcla 

Dependiendo del método que se aplique, se usan las siguientes 
para obtener la viscosidad de la mezcla de fluidos: 

/< ns = ^ A + >g a ' Á \ 


(5.17) 

(5.18) 

(.5.19 

(5 ¿20) 

(.5.21) 

(5.22) 

(5.23) 

(5.24) 

ecuaciones 

(5.25) 


5.6 


FZiij o MíL¿£¿¿ááZco en Tub (lamió 


n — M 
- m 1 » 


“l Ci-h l i 

N 


Donde: 


/U L a* o f o + i“ww 


f = =CL_ 

O B + WOR B 
O w 


f w = - V > 


Donde f Q es la relación aceite -líquido y f es la relación 


agua- 


5.1.8.- Tensión superficial de la mezcla de líquidos. 
Se obtiene con la siguiente expresión : 


<r T - a f + 
l> o o 


W f w 


(5.26) 


C5.27) 


(5.28) 


(5.29) 

•líquido. 


(5.30) 


5.7 


Flujo MuJU¿Í&&¿c.o en TubvtfaÁ 

REFERENCIAS • CAP. 5. 

1. - Beggs, H. D. y Brill, J.P.: "A Stndy of Two Phase Flow in Inclined 

Pipes". Trans. AIME. 1973. 

2. - Mukherjee, H. y Brill, J. : "Liquid Holdup Correlations for Inclined 

Two - Phase Flow". J.P.T. Mayo, 1983. 


5.8 


6. FLOJO MULTIFASICO EN TUBERIAS HORIZONTALES. 
6.1 INTRODUCCICaJ. 


Para flujo horizontal el gradiente de presión debido al cambio de eleva 
cion es igual a cero, por lo que la ecuación (2.9) se reduce á: 


4P— 
A L 


■-AE > 

A L 'f 


+ ( 


AL ac 


0 sea: 


( 


AEL 

AL 


) 

T 


t e v 2 + f A(v^) 

2 d 2 g AL 

c c 


( 6 . 1 ) 


( 6 . 2 ) 


La mayoría de los investigadores han adaptado esta ecuación para apli- 
carla al flujo de dos fases. Para esto suponen que la mezcla gas-Iíqui 
do se puede considerar homogénea en un intervalo pequeño de la tubería. 
Así la Ec. (6.2) se puede escribir como-: 


< -&■ ) 
Ah 'T 


_ÍE_ 


fm 


v 2 

m 


A( V m) 


2 g d 
c 


, 2 g AL 


(6.3) 


En donde f , y v m se refieren a la mezcla y son definidos en - 

forma diferente por los autores de las correlaciones. 

El factor de fricción f , como se indico, depende del número de Reynolds 
esto es, de las fuerzas viscosas y de inercia y de la rugosidad. Para - 
flujo bifásico intervienen además las fuerzas de gravedad e interfacia- 
les. Aunque se ha intentado correlacionar el factor de fricción con gru- 
pos indimensionales que comprenden estas fuerzas, no se ha tenido éxito. 

El enfoque que mas se ha seguido, es determinar el factor de fricción a 
partir de datos experimentales y tratar de correlacionarlo con alguna 
forma del número de Reynolds para dos fases. 

6.2. CALCULO DE LA CAIDA DE PRESION EN TUBERIAS HORIZONTALES. 

Generalmente se considera flujo isotérmico, para el cual las propiedades 
de los fluidos dependen exclusivamente de la presión. En este caso la 
predicción de la caída de presión, consiste en suponer una caída de pre- 
sión A p y aplicar la ecuación (6.3) para determinar el incremento de 
longitud, AL, correspondiente a la Ap supuesta, repitiéndose el pro- 
cedimiento hasta alcanzar la longitud total. Naturalmente que la exacta 
tud de los cálculos aumenta al reducir el incremento de presión supues- 
to? pero también aumenta la cantidad de trabajo requerida para el cálcu- 
lo. Por esta razón se debe establecer un compromiso entre estos dos fac- 
tores, teniendo en mente que el incremento de presión debe ser pequeño 
a presiones bajas, en las que la velocidad varía mas rápidamente con la 
presión, no así a presiones altas en las que la variación es menor; Una 
regla establecida es usar incrementos de presión menores que el 10% del 
valor de la presión media. 
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Cuando no se considera flujo isotérmico, el cálculo del gradiente de - 
presión implica un proceso iterativo, ya que la temperatura es una fun 
cion de la distancia. Entonces, ademas de suponer una A p, se tiene ~ 
que suponer una AL y de ahí determinar la temperatura media de flujo. 

6.2.1. Procedimiento general de cálculo. 

A. continuación se presenta el procedimiento de 'cálculo general para el 
caso de flujo isotérmico. Los pasos 5 a 7 dependen del método que se es 
té aplicando para el cálculo del perfil de presión. “ 

1. - Se inicia con una presión p a la entrada de la tubería. A este - 

punto le corresponde una L 0. 

2 . - Suponga una caída de presión A P y calcule p y p 2 

P * P x “ Ap/2 P 2 = P]L - AP 

3. - Determine las propiedades de los fluidos (R , a , B 0 , Z, B g , ^ o, 

jUg y a las condiciones inedias le escurrimiento . 

4. - Calcule las velocidades superficiales y los gastos de masa de los 

fluidos. Calcule también el colgamiento sin resbalamiento. 

5. - Determine el colgamiento H y la densidad de la mezcla. 

L 

6. - Si las pérdidas de presión por aceleración no se consideran des- 

preciables determine su valor. 

7. - Obtenga el valor del factor de fricción de dos fases. 

8. - Aplicando la ecuación correspondiente determine el valor del gra- 

diente de presión ( Ap/ A L) y con éste la A L correspondiente a 
la A P supuesta. 

9. - Reemplace L por L + A Lj si este valor es menor que la longitud - 

total, haga p¿ - p 2 y regrese al paso (2). Si L es igual o i mayor 
que la longitud total , se termina el cálculo, obteniéndose la pre- 
sión final por interpolación si es necesario. 

Como se indicó este procedimiento se aplica cuando el flujo es isotér- 
mico. Cuando esta condición no es válida, se tienen que incluir los si^ 
guiantes pasos: 

2 1 .- Suponga un incremento de longitud correspondiente a la A p supues_ 
ta y calcule la temperatura media en el incremento. 

8*.- Si la AL C calculada es igual o difiere de la AL S supuesta en 
un valor menor que una tolerancia, continué en (9) . Si no, haga 
AL g = A L , determine la temperatura media en el intervalo y re- 
grese a (3j . 
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6.3.- CORRELACION DE BEE^TIIZZI, TEK Y POETTMANn/ 1 * 

6.3.1. Introducción . 

Las características principales de esta correlación son: 

a) Es independiente del patrón de flujo. 

b) No se consideran las pérdidas de presión por aceleración. 

c) La densidad y el gasto másico de la mezcla están definidos oor las 
ecuaciones : 


fns = + ) 


w m = W L > W g 


d) El factor de fricción para dos fases, f , se obtuvo usando 267 da- 
tos experimentales . Se encontró que los Calores resultantes se po- 
dían correlacionar con una función del numero de Reynolds del lí — 
quido y del numero de Reynolds del gas: 


-r tr^-r 


En donde: 


a * ^/(1 + V ) 

b = 1/exp (0.1 


W g / W L 


Los exponentes a y b se seleccionaron arbitrariamente y para satisfa- 
cer la condición de que la función (6.6) tienda al numero de Reynolds 
del gas cuando la fase líquida tiende a cero, y tienda al número de 
Reynolds del líquido cuando la fase gaseosa tiende a cero. 

La correlación para obtener el factor de fricción se muestra en la 
Fig. 6.1, observándose que es una función de 

6.3.2. Aplicación de la correlación. 


La ecuación del gradiente de presión por fricción, es: 

. 174.158 f„ i tr 

. A£L = tp m 

AL p d 5 

(ns a 
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ooro 
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Fig 6.1 Gráfica del factor de disipación de energía (tomada de Bertuzzi y otros)* 
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El nómero de Reynolds del líquido y del gas se obtienen con las expre- 
siones siguientes. 


ReL 


Reg 


22737 L 


d 


22737 — -S L 

d 


( 6 . 11 ) 

( 6.12 


El factor de fricción puede obtenerse con la Fig. 6.1 o empleando las 
ecuaciones siauientes : ' 2 •' 


ecuaciones siguientes: 

Para : . 0¿pí ¿ 500 


log f tp 

- 1.225 - 0.06561 log 0 - 0.37 * 

(6.13) 

Para: 


0 > 10 000 

■3 ' 

109 f t P 

= 0.49 <P - 0.12616 log 0 - 1.702 

(6.14) 

Para: 



500 6 0 ¿ 10 000 


log f t P 

* P 500 “ 0.6561 y + (1.1056 +1.7723 A F) y 2 



- (0.46214 + 0.90817 AF) y 3 

(6.15) 

donde : 

AF 

F 10000 - F 500 

(6.16) 

F 10000 

= log f fcü (V, 0 « 10 000) 

(6.17) 

F 

500 

- log ' ; * t 0 = 500) 

(6.18) 

y 

= log 0 - 2.699 

(6.19) 


La ecuación (6.15) es un polinomio de interpolación entre los valores 
de f para 0 s 500 y 0 =10 000, con la restricción adicional que la - 
pendiente del polinomio es igual a la de los polinomios (6.13) y (6.14) 
en los extremos de la interpolación. 
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6.3.3. Procedimiento de cálculo. * 

La caída de presión se obtiene aplicando los tres primeros pasos indi 
cados en la sección anterior, y a continuación los siguientes: 

4.- Obtenga w m con las Ecs. (5.16) y (5.20) o con las Ecs. (5.21) , 
(5.22) y (5.23) . Obtenga A con la Ec. (5.2) 

5»- Determine el factor de fricción para dos fases : 

a) Calcule el numero de Reynolds del gas y del líquido aplicando 
las Ecs. (6.11 y (6.12) . 

* b > Con las Ecs. (6.7) y (6.8) calcule a y b. 

c) Determine el valor de 0 de la Ec., (6.6) y con este obtenga el 

valor de f tp de la Fig. 6.2 o con las Ecs. (6.13) a (6.19). 

6.- Obtenga f ns con la Ec. (6.4) 

Aplique la Ec. (6.10) para calcular el gradiente de presión y con 
este obtenga la AL correspondiente al incremento de presión su- 
puesto. 

6.3*4. Ejemplo ’ Aplicando el método dé Bertuzzi, determine el 

gradiente de presión para las siguientes condiciones: 


R 

* 2000 

3 

pies /bl 

Vro 

ti 

o 

00 

q o 

= 5000 

bl/dia a c.s. 

P 

- 1275 lb/pg 2 

d 

= 6 pg. 


T 

= 80 °F 

V 

= 0.7 


]iom 

= 15.18 cp 


Use la correlación de Oistein para obtener las propiedades del acei- 
te. (B Q , Rg X • 

Solución : 

Los pasos (1) y (2) no se aplican 
3. u De la correlación de 0istein: 

B q = 1.1867 

R 460.42 

s 
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De la correlación de Beggs y Rob inson: 

¡ti = fi * = 2.237 cp 

r o L 

Del procedimiento descrito en el inciso 1.4.4 
Z = 0.7786 

í 

De la correlación de Lee: 
u = 0.01396 cp 

g 

De las ecuaciones (1.4.3), (1.4.5) y (1.2.12) 

. 3 s . 3 
pie /pie 

lb/pie 3 

Ib /pie 3 


B = 0.009316 
9 

P = 45.763 




5.7407 


4 - Los valore® de , m y . , se obtienen cpn las Ecs. 

(5.24), y A coilAc. fS.2). 

vv = 4.7648 lbm/seg 
9 

W L = 17.6468 lbm/seg 

w = 22.4116 lbm/seg 
m 

A = 0.3172 


5.a.- Según las Ecs. (6.11) y (6.12) , se tiene 

N = 29891 

ReL 

N = 1 293 548 

Reg 

5.b. - Los valores de «P , a y b, de acuerdo con las Ecs. 
y (6.9) son 

<P * 0.27001 

a = 0.2126 

b = 0.9734 


5.c . - 


0 = (1293548) 0 ’ 2126 ( 29891) ' 9734 
0 = 452 598 


.22) a 


(6.7), (6.8) 
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De la Fig. 6.1 con 0 y « p, se obtiene f = Q.Q05209 

tp 

El mismo resultado se obtiene con la Ec. (6.14) 

6.- De la Ec. (6.4) 

. 3 


“ 18.4359 lbm/pie 3 


7.— Con la Ec. (6.10) se obtiene: 

AP _ (174.158) (0.005209) ( 22. 4116 ) : 
AL ~ c — 

(18. 4359) (6) ' 


-á£L 

Al 


= 0.003179 lb/pg/pie 


6.4.- CORRELACION DE EATON, ANDREWS, KNOWLES Y BROWN (3) 

6.4.1. Introducción. 

Esta correlación se desarrollo' a partir de información sobre condicio- 
nes de flujo obtenida en líneas de 2 y 4 pg de diámetro y de 1700 pies 
de longitud, y una tubería de 17 pg. y lo millas de longitud, en ^nde 
se uso aceite, condensado y agua separadamente como fase líquida y gas 
natural como fase gaseosa. 


Para el cálculo del gradiente de presión se usa la siguiente ecuación: 

(6.19) 


-4£_ _ 

AL 


43.539 f tD 


f 


ns (1 - V 


Donde : 


e k = 


A (y 2 ) + iv A (v ) 2 


9266.1 


w 11 




*^J 

| Ap 


( 6 . 20 ) 


El factor de fricción fue correlacionado a partir de información expe- 
rimental. En la Fig. 6. 2 se muestra la correlación establecida, en la 
que el valor de la abscisa es: 

r -i 0 * 5 

_ 22737 [ W q W m] 

¿^25 ( 6 . 21 ) 


/*g d‘ 
el de la ordenada: 
y = 


w 


o.i 


tp 


( 6 . 22 ) 
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Para obtener las velocidades reales del líquido (v L > y del gas (v ) , es 
preciso conocer el colgamiento del líquido (H^) en cualquier parte de 
la tubería. Esto es necesario solo cuando las perdidas de presión por 
aceleración (Ec. (6.20)) no son despreciables. 

El colgamiento se obtiene con la Ec. (5.1) que proporciona valores _re<a 
les de ésta función. La aplicación de esta ecuación es válida para el 
método de Eaton, ya que durante las pruebas de flujo en las tuberías de 
2 y 4 pg. , se midió el colgamiento de líquido, aislando secciones de 
las tuberías mediante válvulas de cierre rápido. 


6.4.2.- Procedimiento de calculo. 

Los pasos 5 a 8 del procedimiento de cálculo general# para la aplica — 
ción de este método, son los siguientes : 


5. - Cálculo del colgamiento, H . Si las perdidas de presión por acele- 
ración se consideran despreciables, no es necesario determinar e-1 colg£ 

miento. De otro modo H se obtiene con la Ec. (5.1). 

L 

2 2 

6 . - Calcule el valor de A(v ) y A (v ), con las ecuaciones siguien 

tes: L 9 


2 2 2 
A(V L > - V L2 - V LX 


. . 2 . 2 2 

A (v ) = v _ - v 

g g2 gl 


(6.23) 


Con estos valores obtenga E con la Ec . (6.20) 


7.- Cálculo del factor de fricción. - 

a) Determine el valor de x con la Ec. (6.21) y con este obtenga el 
valor de la ordenada, y, de la Fig. 6.2, o con las ecuaciones 
siguientes: 

Para: x ^ 60000 

y s 6677920 X (6.24) 


Para: 

2 3 

60000 < X é 819194 - 39981.7 d + 2838. 8d - 73.26 d 
y = t/100 

t = Antilog (2.37354 - 2.10458 r + 0.5757 r 3 - 0.14189 r 

3 4 

+ S (0.46 - 0.93739 r + 0.45966 r - 0.15975 r ) 

o 3 4 

+ S (0.451 - 0.36293 r - 0.19949 r + 0.12835 r )) 

Donde: r = log (x / 10000) 

S = log d 


(6.25) 

(6.26) 

(6.27) 

(6.28) 
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Para X > 819194 - 39981.7 d + 2838.8 d 2 - 73.26 d 3 

y * (21.525 - 1.5934 d + 0.02278 d 2 + 0.00131 d 3 ) X~°* 49 (6.29) 

b) Calcule el factor de fricción: 


tp 


L 

w 

m 


0.1 


(6.30) 


8.- Aplicando la Ec. (6.19) calcule el valor de (A p/ AL) y con éste el 
valor de AL correspondiente a la Ap supuesta,. 

6.4.3.- Ejemplo.- Calcule la AL correspondiente a un decremento 

Ap = 50 lb/pg^, aplicando el método de Eaton para las mismas 
condiciones del ejemplo anterior f cuando p. .*■ 1300 lb/pg^. Use 
los siguientes valores de R g y B Q . 


( lb/pg 

1300 

1275 

1250 


(pie /bl) 

469.67 

460.42 

451.21 


B 

o 

(bl/bl) 

1.1912 

1.1867 

1.1822 


(1) 

- p i 

- 1300 lb/pg 

h 

= 0 



(2) 


Ap = 50 p 

* 1300 - 

50/2 

- 1275 

lb/pg 2 



P 2 

* 1300 - 

50 = 

1250 

lb/pg 2 

(3) 

De 

correlaciones se 

obtiene : 





P 

2 



B 

g 

a 


1300 

0.7749 

2.199 


0.009093 

10.68 


1275 

0.7786 

2.237 


0.009316 

10.87 


1250 

0.7823 

2.276 


0.009548 

11.06 


(4).- A las presiones de interés se obtiene: 


1300 

1275 

1250 


45.664 

45.763 

45.861 


sL 

1.9669 

1.9614 

1.9539 


sg 

4.1011 

4.2271 


6.11 


4.3583 
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y a p ~ 1275 lb/pg , se obtiene: 
f - 5.7407 

v g 

w « 4.7648 

g 

w * 17.647 

L 


- 0.01396 cp 


w * 22.412 
m 


0.3172 


(5) . - Cálculo de H. . y H__ 
L1 L ¿ 


N l = 0.15726 [l / a 3 j0.25 

N l1 = (0.15726) (2.199) (1/(45.664) (10.68) 3 ) °* 25 


N l1 = 0.022518 


N Lv -M 38 * sL <Pl /CT) 


0.25 


N Lv i - (1.938) (1.9669) (45.664/10.68) 


0.25 


N t „ * 5.4813 
Lvl 


gv 


1 * 938 V sg ( Cl / a > 


0.25 


N = (1.938) (4.1011) (45.664/10.68) 


0 . 25 


N _ = 11.429 
gvl 


La ecuación de Mukherjee, oara 9 = 0, es: 

2 N« 

H L " SXP < (C 1 + C 4 \ 

N 

Lv 


H - exp ( (-.38011 + (2.34323) (0.022518) 2 ) - ( — - 4 — - 

(5.4813) ü * 


H = 0.4777 

L1 


47569 

28866 
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De la misma manera se obtiene H. 


L2 


N L2 * °- 022678 

N t 0 = 5.4035 

Lv2 

N - » 12.0529 
gv2 

«L2 = 0.4673 


(6) .- Cálculo de A (v 2 ) , A(v 2 ) y E v 
L g K 


sLl 1.9669 


“ ' H u ‘ 

_ 1-9539 

V L2 0.4673 


o f 4 7 7 7 ~ ^ 4.1174 pies/seg 

= 4.1813 pies/seg 


sgl 


4.1011 

v gl (1 - H l1 “ 0.5223 


«7.8520 pies/seg 


4 ICQI 

v - « 8.1815 pies/seg 

g2 0.5327 


A(v ) 2 = 4.1813 2 - 4.1174 2 « 0.5303 

L 


A (v j 2 = 8.1815 2 - 7.8520 2 « 5.2830 

g 


Sustituyendo valores en la Ec. (6.20) 

(17.647) (0.5303) + (4.7648) (5.2830) 


E„ 


9266. 


17.647 

4.7648 

45.763 

5.7407 . 


(50) 


E r « 0.00006131 


(7) . - Cálculo del factor de fricción. 


(7.a).- Determinación del valor de la ordenada. Sustituyendo valores 
en la Ec. (6.21) : 


X = 


(22737) (4.7648) (22.412) 


0.5 


(0.01396) (6) 


2.25 


« 298724 
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De la Fig. 6.2 


v c -^-) 0 - 1 

tp V 


0.016 

T7¡6j7W ri 

22.412 


fns = Á ?L +(1 - A) fe, 


ns = (0.3172) (45.763) + (0.6828) (5.7407) 


P - 18. 436 
\ ns 


Sustituyendo valores en la Ec. (6.19) 


A P ^ (43.539) (22.412) (0.01639) 

AL (6) 5 (18.436) (1 - 0.00006131) 


0.0025 lb/pg 2 /pie 


(9) .- AL = Ap/dp/AL * 50/0.0025 = 20 000 pies 


6.5.- CORRELACION DE BEGGS Y BRILL 


6.5.1.- Introducción. ■ ' 

Esta correlación fue desarrollada a partir de datos de flujo obtenidos 
en tuberías de acrílico de 1 y 1 1/2 pg y de 90 pies de longitud, las 
cuales se podían inclinar a cualquier ángulo. Los fluidos utilizados 
fueron aire y agua. 

La siguiente expresión se emplea para calcular el gradiente de presión: 


-4? — = 43.539 f tp W m 


7.2557 fm W g 

fn S P d4 K 


Definiendo: 


7.2557 W m 

6is Pd4 fg 


La Ec. (6.24) se transforma a: 

AP . 43.539 f t p w m 

fns * 5 (1 -V 
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El factor de fricción se calcula de: 
f . 


f * ( — \ f 
tp v f ' n 

n 


Donde f fl es el factor de fricción del diagrama de Moody para tuberías 
lisas. Beggs da la siguiente expresión para su obtención: 


f n = _ 2 [ t 


5223 log N - 3.8215 
Re 


En donde: 


124 d v m (°n 
Pns 


“ns = “l 1 + >‘ g (1 ' A > 


El factor de fricción para dos fases normalizado (f / f ) se encon^- 
tró ser función, del colgamiento de líquido H y de p A n , y se 
puede calcular de: , 


- 0.0523 + 3, 182 lny - 0.8725 (lny ) 2 + 0 . 01853 (lny) 4 (6.39) 


Las tuberías transparentes permitieron observar y clasificar tres pa- 
trones de flujo* De acuerdo a esta clasificación se desarrolló el ma- 
pa de patrones de flujo, mostrado en la Fig. (6.3) , en función de X 
y del numero de Froude. El patrón de flujo para cualquier condición 
de flujo se puede determinar de este mapa o de la tabla siguiente: 


Patrón de Flujo 


Condiciones 


Segregado 


x < 0.01 y N pR * L x Ó A ¿ O.oi y N pR ^ L,, 


Transición 


a >0.01 y L 2 C N fr ^ L 3 
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Número de Froude 



0.0001 0.001 

Colgamiento sin resbalamiento, X 

Fig 6.3 Mapa de patrones de flujo de Beggs {4) 


6.16 



Flujo UvúU¿Ífa¿CLO m TubeAúu Ho/Uzonlaleó 


Intermitente 0.01 ^ A iz 0.4 y L_ < N ^ L, 6 

3 FR 1 


A> 0.4 y L, < N 

3 FR 


Distribuido 


A < 0.4 y N pR ^ L, 6 A > 0.4 y N M > L 


FR 


En donde: 


N * 12538 
FR 


2 .5 




(6.41) 


0.302 


L x * 316 A 


L 2 = 0.0009552 A 


-2.4684 


L 3 * 0.1 A 


-1.4516 


L . ■ 0.5 A 

4 


-6.738 


El colgamiénto se obtiene con la Ec. (5.2) 

6.5.2. Procedimiento de Cálculo . - 

5. ~ Calcule el colgamiento con la Ec. (5.1) y f con la (5.9) 

6. - Determine el valor de E^ con la Ec. (6.32) . Sí se consideran - 

despreciables las pérdidas de presión por aceleración, haga 

E =0. 

K 

7. - Con las Ecs. (6.34) a (6.37) determine (f^/f^) y f ft 


8.- Calcule f , de f * (f / f ) f 
tp tp tp n n 


9.- Obtenga ( APAL) aplicando la Ec. (6.33) y con este valor la 
AL correspondiente a la. AP supuesta. 


Ejemplo .- Aplicando el método de Beggs y Brill, calcu- 

le el valor de la AL correspondiente a un decremento de presién 
Ap » 50 lb/pg 2 , para las mismas condiciones del ejemplo anterior. 

(1) P x » 1300 lb/pg 2 

(2) p ■ p 1 - Ap/2 = 1300 - 50/2 * 1275 lb/pg 2 
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(3) Del ejemplo anterior se tiene: 

R s = 460.42 pies 3 /bl 

B = 1.1867 bl/bl 

o 

z = 0.7786 

B g = 0.009316 pie 3 /pie 3 

fi = 2.237 cp 

i Xí 

= 0.01396 cp 

Cf = ' 10.87 dinas/cm. 

(4) f> L = 45.763 lb m /pie 3 

O = 5.7407 Ib /pie 3 

Vg m 

w = 4.7648 Ib /seg 

g m 

w = 22.412 Ib /seg 

m m 

v = v + v =1.9614 + 4.2271 = 6.1885 pies/seg 
m sL sg 

A = 0.3172 

P 18.436 Ib /pie 3 

v ns m 

(5) 

N = 0.15726 (1/ P T (T 3 ) 0.25 

L 1j Jj 

0.25 

N = (0.15726) (2.237) (1 / (45.763) (10. 87) 3 ) 

L 

N = 0.022593 

L 

N = 1.938 V T C f T / <J}°* 25 

Lv sL v L 

0.25 

N = (1.938) (1.9614) 045.763/10.87) 

Lv 

N = 5.4441 

Lv 

N = 1.938 v ( p./CÍ) 0 - 25 

gv sg v » L 
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N = (1.938) (4.2271) (45.763 / 10.87) 
gv 


N = . 11.7346 
gv 

' -r . 7346) 0,47569 

H = exp (- . 38011 4 (2.34323) (0.022593) 2 ) r- ‘ A' 

L V* ( 5.4441) 


H = 0.4725 

Xi 


e m = í*l b l 


(í - h l j e 


p * (45.763) (0.4725) + (0.5275) (5.7407) 
v m 

P * 24.65 Ib /pie 3 
V m ni 

(6) Con la Ec. (6.32) se obtiene: 

_ (7.2557) (24.65) (22.412) (4.7648) 

E K “ (18.436) (1275) (6 4 ) (5.7407) 


E k = 0.0001092 


(i - X ) 


= (2.237) (0.3172) + (0.01396) (0.6028) 


u . * 0.718 cp 

ns 


(124) (6) (6.1885) (18.436) 


De la Ec. (6.35) se obtiene; 


= 1.4202 


El valor de S, según la Ec. (6.39) es: 


S *■ 0.36662 
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tp s 0.3666 

= e = e * 1.4428 


f tp = f n e =(0.017393) (1.4428) =0.025095 


( 8 ) 


AP _ (43.539) (0.025095) (22.412) ; 


AL 


(18.436) (6) ( 1 -0.0001092) 


_A£ 


Jl = 0.003829 lb/pg /pie 


(9) 


AL = ( Ap/ AP / AL) = 50/0.00382 = 13 058 


A continuación se presentan los valores de Ap/ AL obtenidos median 
te la aplicación de los tres métodos explicados, para el mismo proble 
ma: 


Método 

Bertuzzi, Tek y Poettmann 
Eaton, Andrews, Knowles y Brown 
Beggs y Brill 


AP / AL 
0.003179 
0.0025 
0.003829 
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7,~ FLUJO MULTIFAS ICO VERTICAL 
7.1.- Introducción. 

Es evidente^ la importancia de la evaluación de las caídas de presión - 
en la tubería vertical, ya que la mayor proporción de la presión dis- 
para llevar los fluidos del yacimiento hasta los separadores 
se consume en dicha tubería. 

En^la siguiente tabla se muestran las proporciones en que cae la pre- 
sión en el yacimiento , la tubería vertical y la línea de descarga del 
pozo, para cuatro valores de índice de productividad y gastos de acei- 
te. 

DISTRIBUCION DE PERDIDA DE PRESION EN FLUJO DE ACEITE (1) 

CARACTERISTICAS 
DEL POZO 


D = 10000 pies 
P ws i 3000 UVpg 2 

R = 750 pies 2 /bl 

9 TP “ 3 V2 pg 

P th - 100 lb/pg 2 

Dada la magnitud de las pérdidas de presión en las tuberías de produc- 
ción se hace indispensable su evaluación precisa, a fin de optimizar el 
sistema de producción de los pozos. 

La determinación de las distribuciones de presión en las tuberías de - 
producción permites 

a) diseñar las tuberías de producción y líneas de descarga 

b) Obtener el punto óptimo de inyección de gas en el bombeo neumático. 

c) Proyectar aparejos de producción artificial, (neumático, mecánico, 
eléctrico) . 

d) Obtener p wf sin necesidad de intervenciones en los pozos. 

7.1.1.- Flujo multifásico en tuberías verticales. 

Cuando fluyen simultáneamente aceite y gas a través de una tubería ver- 
tical a medida que se incrementa la proporción de gas en el flujo, las 
caídas de presión tienden a disminuir, hasta alcanzar un mínimo. A con 
tinuación los aumentos en la cantidad de gas provocaron incrementos en 


J 

q o 

% DE PERDIDA I 

(bl/dla/lb/pg 2 

(bl/día) 

1 

Yacto. 

T.V. 

L.D. 

2.5 

2700 

36 

57 

7 

5.0 

3700 

25 

68 

. ■ 

7 

10.0 

4500 

15 

78 

7 

15.0 

4800 

11 

82 

7 
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las perdidas de presión,. Fig. 7.1 (a) . Este fenómeno se explica de la - 
manera siguiente: 

a) Para volúmenes pequeños de gas plevalece la carga de líquido, misma 
que va reduciéndose al aumentar dicho gas; ya que la densidad de la mez^ 
cía gas-líquido continuamente disminuye. 

b) Después de que el volumen de gas alcanza cierta proporción, las pér- 
didas por fricción debidas al flujo del propio gas aumentan notablemen- 
te, compensando y sobre pasando la disminución en la carga hidrostática. 
El efecto resultante es el aumento en las caídas de presión. 

Asimismo , manteniendo los gastos de líquido y gas y variando el diáme- 
tro del conducto, se ha observado un comportamiento similar al descri- 
to; conforme se aumenta el diámetro, primero disminuyen las pérdidas de 
presión hasta un mínimo, y luego aumentan indefinidamente. Fig. 7.1 (b) . 
Eh este caso el mecanismo que prevalece después del mínimo es el resba- 
lamiento entre las fases; ésto es, que el gas viaje a una velocidad ma- 
yor que la del líquido, lo que implica un retraso de éste respecto al 
gas, resultando en mayor carga hidrostática. Para diámetros grandes de 
tuberías, la velocidad del líquido es baja y el retraso entre las fases 
es notable. Adl disminuir el diámetro aumenta la velocidad del líquido; 
y aunque la del gas también aumenta, lo hace en menor proporción por su 
compresibilidad , el resultado es una reducción en la caiga hidrostática. 

El efecto del resbalamiento se visualiza más fácilmente observando lo 
que ocurre en un tanque cilindrico lleno de líquido, al que se le está 
burbujeando gas en el fondo. Evidentemente las burbujas de gas de segre_ 
garán del líquido liberándose en la superficie. Supóngase ahora que se 
va reduciendo él diámetro del tanque. Se alcanzará un diámetro en que 
el gas ya no resbale y empiece a arrastrar parte del líquido existente. 

¡ . 

Por otra parte si se mantiene fijo el gasto de gas en un conducto verti- 
cal y se varía el volumen de líquido se tendrá por efecto, del resbala- 
miento el siguiente comportamiento: 

a) Para bajos gastos de líquido el resbalamiento será grande y la dife- 
rencia de presiones entre dos puntos del conducto, se deberá principal- 
mente a la carga de líquido. 

b) Al aumentar el gasto de líquido tenderá a disminuir el resbalamiento, 
lo que se traducirá en la disminución en la carga de líquido y una re- 
ducción en las pérdidas de presión. 

c) Para gastos grandes de líquido las pérdidas por fricción compensarán 
la reducción de la carga hidrostática, incrementándose las caídas de pre 
sión. 

Los tres casos anteriores se ilustran en la Fig. 7.1 (c) 

Debido al resbalamiento no es posible calcular la relación gas-líquido 
a condiciones de flujo a partir de las condiciones de entrada a la tube_ 
ría. Dicha relación se obtiene a través de correlaciones, que han sido 
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Fig 7.1 Relación entre la caída de presión y q R y <¿> , en 
tuberías verticales 



Fig 7.2 Mapa de patrones de flujo de Duns y Ros con modi- 
ficaciones 13 
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desarrolladas experimentalmente en base a la distribución de las fases 
en la corriente. Por lo general, de las correlaciones citadas se deter- 
mina el valor del colgamiento "HOLD-UP" y de éste, fácilmente pueden ob 
tenerse los volúmenes de gas y líquido en el flujo. 

7.2. - Patrones de flujo. 

( 2 ) > 

Ros N.C.J. identifico seis patrones típicos de flujo multifásico en 
tuberías verticales que denomino: burbuja, tapón, bache, espuma, tran- 
sición y niebla; sin embargo en la mayoría de las correlaciones estable 
cidas no se consideran los regímenes de flujo tapón y flujo de espuma. 

Ros observo las siguientes condiciones de flujo: 

a) Para bajos gastos de gas prevalece el flujo de burbuja, la fase lí- 
quida es continua y el gas esta disperso en burbujas pequeñas. 

b) A mayores gastos de gas, pero gastos bajos de líquido, conforme au- 
menta el gas, el numero y tamaño de las burbujas también aumenta, 
tomando forma de bala (flujo tapón). A continuación estas burbujas 
coa leseen formando baches que contienen principalmente gas y que al- 
ternan con baches de líqudio (flujo de bache). 

c) Para v sg » 50 pies/seg y v s l < 1.25 pies/seg el flujo cambia de ta 
p5n a niebla. 

d) Cuando v sL alcanza valores superiores a 5.25 pies/seg ya' no es fácil 
distinguir los varios patrones de flujo. 

e) Para valores bajos de v sg y v sL se presenta el fenómeno conocido co- 
mo cabeceo en el que el flujo varía cíclicamente en pocos segundos. 

El flujo es inestable y los gradientes de presión son muy variables 
y difíciles de predecir. 

El patrón de flujo existente dentro de la tubería vertical se obtiene, 
generalmente, en función de los números adimensionales de velocidad del 
gas y del líquido, o en función de las velocidades superficiales. Los 
autores que optaron por estos enfoques desarrollaron mapas de patrones 
de flujo en los que se delimitan las regiones en que ocurren los regí- 
menes considerados. En las Figs. 7.2 y 7.3 se presentan dos ejemplos de 
mapas de patrones de flujo. 

7.3. - Enfoques en el desarrollo de las correlaciones. 

Las diversas correlaciones existentes para el cálculo de distribuciones 
de presión con flujo multifásico, pueden clasificarse en, tres tipos bien 
definidos: 

Tipo I.- No se considera resbalamiento entre las fases. La densidad de 
la mezcla se obtiene en función de las propiedades de los fluidos, corre^ 
gidas por presión y temperatura. Las perdidas por fricción y los efec- 
tos del colgamiento se expresan por medio de un factor de fricción co- 
rrelacionado empíricamente . No se distinguen patrones de flujo. Dentro 
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Fig 7.3 Mapas de patrones de flujo de Chierici , Ciucci y 
Sclocchi ,, 
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de este tipo están incluidos los métodos de Poettmann y Carpenter 
cher y Brown' 4 ' y Baxendell y Thomas. 


í 3) 


Fan 


Tipo II.- Se toma en cuenta el resbalamiento éntrelas fases. La densi- 
dad de la mezcla se calcula utilizando el concepto de colgamiento. El 
factor de fricción se correlaciona con las propiedades combinadas del 
gas y líquido. No se distinguen regímenes de flujo. El método de Ha- 
gedorn y Brown' 6 ' cae dentro de este tipo de correlaciones. 

Tipo III.- Se considera resbalamiento entre las fases. La densidad de 
la mezcla se determina mediante el colgamiento. El factor de fricción 
se correlaciona con las propiedades del fluido en la continua. Se 

distinguen diferentes patrones de flujo. Las principales correlaciones 
que caen dentro de esta clasificación son: Duns y Rostí), orkiszewskit 7 ) , 
Azizt Beggs y Brill^ 9 ), Chierici f 10 \ , Gould y Tek f tH), etc. 


7.4.- METODO DE POETTÍÍANN Y. CARPENTER 
7.4.1. Introducción . 


(3) 


Poettmann y Carpenter publicaron en 1952 un procedimiento analítico para 
determinar las caídas de presión en tubería verticales con flujo multi- 
fásico. Su ecuación principal la desarrollaron a partir de un balance 
de energía entre dos puntos dentro de la tubería de producción. Esta 
ecuación es : 


dh 


1 

144 


_ + 
ns 


£ tp. (q o M) ^ 


(7.1) 


2.979 x 10" 


P d 
ns 


El factor de fricción fue determinado a partir de datos medidos de pre- 
siones de fondo en 49 pozos fluyentes y del bombeo neumático y aplican- 


do la ecuación anterior. Los valores de 
naron con el numerador del numero de 
dades prácticas queda: 


así obtenidos se corre lacio- 
que expresado en las uni 


d v p = 1 . 77 x "lÓ 

ns 


“o M 


(7.2) 


■ (13) 

La correlaciórí así obtenida fue' extendida por Baxendel y Thomas , pa 

ra ser aplicable a pozos con altos gastos y flujo por’ el espacio anular. 
En la Fig. 7.4 se muestran los resultados de ambas correlaciones. 

La siguiente ecuación puede emplearse para obtener el valor de f 


f. = 5.415 x 10" 3 
tp 

+ 3.5843 x 10 


5.723 x 10’ 4 a +1.848 x lo” 4 a 2 


-6 


(7.3) 
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Fig 7.4 Correlaciones del factor de friccio'n 1 



Fig 7.5 Correlacio'n para el factor de friccio'n, Fancher 
y Brown (4) 
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Donde : 


d x 10 

q o H 


(7.4) 


7.4.2. Procedimiento de cálculo. (Fig. 7.6) 

1. - A partir de una p y h, dadas (condiciones en la cabeza o en el fon 

do del pozo) , fijar una A p y obtener 

P 2 = P x + YP S P x + ( dp/2) 

2. - Calcular para las condiciones medias del intervalo p, T (esta tempe 

ratura generalmente estimada) , los valores de Z, B 0 , R s y 3 g . 


3.- Calcular a p con la Ec. (5.18) 


4.- Determinar el valor de d v p y obtener f de la Fig. 7.4 
o con las Ecs. (7.3) y (7.4)^ S P 


5. - Aplicando la Ec. (7 . 1) , cuantificar á h 


6.- Repetir el procedimiento hasta completar la profundidad total del 
■ . pozo. . 


(4) 


Siguiendo el mismo enfoque de Poettmann y Carpenter , Fancher y Brown 
establecieron una correlación para determinar el factor de fricción en 
función de la relación gas líquido producida. Fig. 7.5. 


7.4.3.- Ejemplo.- Se desea calcular el gradiente de presión en un pozo 
con flujo vertical bajo las siguientes condiciones : 


q * 943.5 bl/día 

° 3 

R - 1.122 pies /bl 


d = 1.995 pg 
p = 960 lb/pg 2 


T - 153°F 


y =0.85 
r ro 


y = 0.65 

9 


2.- Obtener las propiedades de los fluidos a las condiciones de presión 
y temperatura : 

Con la correlación de 0istein: 


De las otras correlaciones, se obtiene: 
Z = 0.903 


R = 212.19 
s 


B = 0.01628 
9 
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FIG. 7.6 DIAGRAMA DE FLUJO SIMPLIFICADO, PARA FLUJO MULTIFASICO 
EN TUBERIAS VERTICALES. 

METODO DE POETTMAN Y CARPENTER. 
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3.- De la Ec. (5.18) : 


. 3 


p =16.825 lbm/pie 
r ns 


4. - De la Ec. (7.4) 


Con la Ec. (7.3) 


a = 5.98 


f = 0.00938 
tp 


5.- Aplicando la Ec. (7.1) , se obtiene: 


= 0.1625 lb/pg 2 /pie 


7.4.4.- Obtención del gasto optimo. 

Según se indico, la caída de presión en una tubería de producción, es mí 
nima para un gasto determinado. Fig. 7.1 (c) . Este gasto ha sido de fin i”* 
do como el gasto optimo. Para gastos mayores la caída de presión aumenta 
por el efecto de la fricción. Para gastos menores el colgamiento del lí- 
quido origina el incremento en la caída de presión. 

Para gastos bajos la acumulación de líquido provoca un incremento consi- 
derable en el peso de la columna de fluidos. Este aumento en la carga de 
fluidos, reduce la velocidad del flujo, lo que a su vez causa un mayor 
resbalamiento. El resultado de esta secuencia es la precipitación del 
flujo en un estado inestable, que produce rápidamente la suspención de 
dicho flujo. 

El gasto optimo puede obtenerse aplicando la correlación de Poettmann - 
Carpenter - Baxendel - Thomas , de la manera siguiente: 


Sustituyendo (7.4) y (7.3) en (7.1), se tiene : 


;áJL. = 
dh 


144 


p + 

ns 


( q o M ) 2 

2.979 x 10 5 P~ 
ns 


[>• 


415 x 10 


-3 


5.723 x Kf 4 ( -AJLJO 


+ 3.5843 x 10 


-6 


q Q M 


d x 10 


- ) + 1.848 x 10~ 4 ( 


q M 


d x 10 2 

q M ■ 
o 


(7.5) 


Derivando con respecto a q Q M e igualando a cero: 


d ( dp/ dh ) 
d ( q M ) 


144 x 2.979 x 10^ d 5 p 


-v~3 


2 x 5.415 x 10 ’ q M 
OP 
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- (5.723 x 10’ 4 ) (d x 10 6 ) - J 1, 5843 x 1Q~ ) (d x 10 V 


( q__ m* 

°p 


= O (7.6) 


Simplificando ; 


íq op M) 3 5 * 2844 x 1q4 d (q Qp M) 2 “ 3.3096029 x 10 14 d 3 = O (7. 7) 


Resolviendo la ecuación de 3er. grado anterior, para d = 1, 2, 3, 4 y 6 


pg. , se obtiene : 



Para: d = 1 

pg- 

q M = 
o 

91970 lbm/dia 

d = 2 

pg. 

li 

183942 

d = 3 

pg* 

q M = 
o 

275913 

d = 4 

pg. 

q o M = 

367883 

id 

ll 

p 

pg- 

q Q M = 

551825 

Los valores anteriores 

constituyen una recta, 

cuya ecuación es 


q M = 91970 d (7.8) 

op ' 

Sustituyendo M según la Ec. (5.16) y despejando q , se tiene: 

op 

q__ = 91970 d / (350.5 V + 0.0764 R y + 350.5 V WOR) ( 7 . 9 ! 
op ro ' g w 


La ecuación anterior porporciona un medio sencillo y útil para diagnosti 
car las condiciones de flujo en los pozos y para seleccionar las tube- 
rías de producción. 

Ejemplo.- Se desea diagnosticar las condiciones de flujo actuales y se- 
leccionar la tubería de producción óptima para los pozos productores de un 
campo, de los que se ha obtenido la siguiente información: y - 0.80; 

R = 2000; WOR =0; y =0.75. La producción promedio de los r ° pozos es de 
250 blo/dia, f luyendo 9 pqr una T.P. con d = 2.441 pg. 

a) Diagnóstico.- Sustituyendo datos en la Ec. (7.9) se obtiene el gasto 
óptimo para las condiciones de flujo actuales: 

q op = {91970) (2-441) / ((350.5) (0.80)+(0.0764) (2000) (0.75)) 

g = 568 blo/dia 

"op 


Como el gasto real es mucho menor que el óptimo, se infiere que el flujo 
es ineficiente, debido a un colgamiento excesivo que incrementa el gra- 
diente de presión. 

b) Selección de la tubería óptima.- La misma Ec. (7.9) puede emplearse 
para obtener el diámetro óptimo, suponiendo un gasto de 250 blo/dia 
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d .= q o M / 91970 = 1.073 pg 

Aunque el método es aproximado, permite en el caso analizado, con- 
cluir la conveniencia dé usar tuberías de menor diámetro, a fin de 
prolongar la vida fluyente de los pozos. 

7 . 5 . - METODO DE BÉÓGS Y - BRILL . 

7.5.1.- Introducción . 

Estos autores establecieron, a partir de pruebas de laboratorio , una co- 
rrelación para calcular las distribuciones de presión en tuberías con 
flujo multi fásico. El procedimiento es aplicable a flujo horizontal, in- 
clinado y vertical. / , ' - — - o 

La experimentación fue realizada en tubos transparentes de acrílico de 
90 pies de longitud. Estos tubos fueron dotados dé un mecanismo que per 
mitio variar su posición desde la horizontal a la vertical; además se 
les incluyeron dispositivos para medir gastos, caídas de presión, ángu- 
los de inclinación' y* el colgamiento . Los fluidos utilizados fueron aire 
y agua. Las pruebas consistieron en medir el colgamiento y la diferen- 
cial de presión en una sección de la tubería. En total se tomaron 584 
datos de colgamiento y caída de presión. 

Dichas pruebas se realizaron para los siguientes rangos: 1) Gasto de gas : 
0 a 300 m pies 3 /dí a , 2 ) Gasto de líquido : 0 a 1030 bl/día; 3) Presión 
media del sistema: 35 a 95 lb/pg . 2 abs. ; 4) Diámetro de 1 y 1.5 pg. ? 

5) Colgamiento de 0 a .87 y 6) Angulo de inclinación: -90 a +90°. 

No obstante que el método fue desarrollado dentro de. rangos limitados?*, 
en trabajos posteriores se ha comprobado que permite predecir con bas- 
tante exactitud las caídas de presión , en tuberías verticales con flujo 
simultáneo de aceite, gas y agua. 

Los mismos autores integraron un banco de datos con información de más 
de 700 pozos, mismo con el que evaluaron los resultados que se obtienen 
calculando las caídas de presión con diversos métodos, incluido el de 
ellos. 

A pesar.de que un método se ajusto mejor a los datos considerados, no es 
posible seleccionar un método para aplicarlo en la generalidad de los ca 
sos. Ante esta situación los datos de agruparon por rangos de : diámetro” 
de tubería, densidad API, relación gas-líquido, relación agua-aceite, ve 
locidad superficial y se establecieron comparaciones entre los resulta-” 
dos calculados y los medidos? pero tampoco fue factible determinar cual 
método es el más adecuado para ciertas propiedades. 

La ecuación general establecida es: 


dp 

_ 

-2_ 

g c 

, i 

con Q a 4. — 

tp_ 

_ 2 
p V 

r ns m 

OCll ü y ; 

m • 


5.362 d 

dz 

144 



~ 





i _ l£SL 

—BL 

- 




g c . . Vi 

? 
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Cuando H 1, se reduce a la ecuación para la fase líquida. 

L 

Cuando H 0, se reduce a la ecuación para la fase gaseosa. 

L 

Cuando 0 = 0, se reduce a la ecuación para flujo horizontal. (Ec. 6.33) 

Cuando 9 * - 90°, se reduce a la ecuación para flujo vertical. Para e£ 
te caso la ecuación anterior es: v 



En esta ecuación se identifican los gradientes por densidad y por fric 
ción, así como el termino de aceleración, o sea: 


dh 


v dh ’ e ■> dh /f 


(7.12) 


7.4.2. Procedimiento de cálculo. 

Después de los tres primeros pasos, se continua con: 

4. - Calcular v__, v ' v y A , con las Ecs. (5.3), (5.4), (5.5) y (5.6) 

SL, m 

5. - Obtener el colgamiento, H , con la Ec. (5,1) 

ÍJ 

6. - Calcular p con la Ec. (5.9) y p con la Ec. (5.18) 

m r ns 

7. - Determinar u con la Ec. (6.37) y N con la (6.36) 

• ns «e 

8. - Calcular f con la Ec. (6.35) , y f con las Ecs. (6.38), (6.39) y 

(6.40). P 

9. - Obtener el término de aceleración E^ 

10. - Calcular el gradiente de presión total con la Ec. (7.13) y el valor 

de Ah. 

11. - Repetir el procedimiento hasta completar la profundidad total del 

pozo. 
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Ejemplo.- Calcule el gradiente de presión para las mismas condiciones 
de flujo del ejemplo anterior. Considere despreciable el efecto por - 
aceleración. El valor de a f es de 10 cp. 

3. Del ejemplo anterior y de otras correlaciones, se tiene; 


* 1.1002 

* 212.19 

>■' 3.075 
. - 0.06848 


- 0.903 
» 0.01628 

= 13.21 
» 49.916 
= 3.0495 


4. v sg * 7.4509 


=10.5572 


A = 0.2942 


5. Obtención de H . (Ec. (5.1) ) 
L 

N l = 0.026256 


N T ■ ,»■ 8.3933 

LV 


N = 20.132 

gv 


Sustituyendo: 


> ( - 0 . • 


H L * exp ! ( ~ 0*38011 + 0.12988 - 0.11979 + (2. 34323) (0.026256 ) Á ) 


20.132 


0.47569 


8.3933 


0.28866 


= 0.43539 


6 . 


P ■ (49.916) (0.43539) + (3.0495) (0.5646) = 23.455 


P__~ 16.825 

ns 


7. Sustituyendo datos en las Ecs. (6.37) y (6.36) 

^ns = ’.K*- + *g (1 


(3.075) (0.2942) + (0.06848) (0.7058) = 


7.14 


0.9530 
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124 d V m P T 


(124) (1.995) (10.557) (16.825) 
(0.9530) 


N_ - 46 107 

Re 


8. De la Ec. (6.35) , se obtiene: 


f = 0.021293 


El valor de y, según la Ec. (6.40) es: 


= 1.55198 


(0 . 43539) 


De la Ec. (6.39) se obtiene: 


S = 0.37298 


De la Ec. (6.38) : 


f * 0.021293 e 
tp 


= 0.030919 


10. Sustituyendo valores en la Ec. C7.ll 


dp 1 . . (0.030919) (16.825) (10.5572) 

5T = 144 ( 23 ' 455 + (5.362) (1.995) 1 


= 0.20052 lb/pg 2 /pi 
dn 


7.6.- FLUJO EN CONDUCTOS ANULARES. 


Los métodos descritos anteriormente pueden aplicarse al flujo multi fá- 
sico por el espacio anular. Para esto basta sustituir el diámetro de 
la tubería por el diámetro hidráulico (d^) y considerar el área anular 
(Ap) en lugar de la de la tubería, o sea: 


4 área de la sección de tubería 
n ” perímetro mojado 

» < d ci - 4 > 
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Para un conducto anular: 



(d 2 . - d 2 
ci te 


* (d . + d ) 
ci te 



Para flujo por el espacio anular se ha considerado que el concepto de - 

diámetro hidráulico es válido cuando d /d -<0.3 

te ci 

( 12 ) 

Cornish obtuvo resultados muy satisfactorios, al aplicar el método 
de Poettmann - Baxendel, usando como rugosidad absoluta el termino si- 
guiente : 


r d . 

C1 



d te I 

d . + d^_ 

L ci te 


te 

d . + • d, 
ci te 


Donde: 

'f C £-~ rugosidad de la superficie interior de la T.R. , pg. 

C1 ( *ci = 0.0018) 

€ fce .- rugosidad de la superficie exterior de la T.P., pg. 

( 6 te « 0.006) 

7.7 SELECCION Y AJUSTE DEL METODO DE CALCULO ‘ 

Es evidente que para calcular las caídas de presión por las tuberías, - 
se seleccionará el método que permita reproducir, con mayor aproxima- 
ción, los valores medidos en el campo, de dichas caídas de presión. 

Sin embargo generalmente es necesario realizar algunos ajustes para lo- 
grar mayor precisión en los resultados . Los valores de las propiedades 
de los fluidos deben corregirse, según se indicó en el capítulo 1. Si 
se observa una desviación sistemática en cierto sentido, entre los va- 
lores calculados y los medidos, entonces el ajuste es sencillo. Basta- 
rá en este caso incluir un factor de eficiencia o modificar el valor 
de la rugosidad de la tubería, hasta lograr la mejor aproximación. 

7.8 FLUJO EN BACHES 

En las tuberías verticales submarinas ("risers") , que descargan en las 
plataformas de producción, frecuentemente el flujo se presenta en forma 
de baches. La presencia de éstos ocasiona que los separadores de gas y 
aceite y las bombas operen de manera ineficiente. Este tipo de flujo pue 
de eliminarse reduciendo el diámetro del "riser M : sin embargo se origi- 
na una contrapresión que provoca la disminución en la capacidad de pro- 
ducción del sistema. 

Schmith^^ investigó la influencia del flujo, tanto en la tubería hori- 
zontal de llegada como en el "riser" , en la formación de los baches. 
Observó que el flujo de baches en los "risers” ocurría cuando en la tu- 
bería horizontal de llegada el flujo era estratificado o en baches. 
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Respecto a los patrones de flujo, determino que la correlación de Beggs 
y Brill es aceptable para las tuberías horizontales. Los mapas propues- 
tos por Duns y Ros^ y Aziz^ 8 ^, fueron los mejores para predecir el 
tipo de flujo en los "risers". 

Para evitar los efectos indeseables del flujo en baches se ha sugerido 
el empleo de un "riser" múltiple, constituido por varias tuberías con 
un área equivalente a la de la tubería sola. También sp ha propuesto 
reinyectar el gas obtenido del separador, para romper el bache y cam- 
biar el patrón de flujo a tipo espuma. 
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8,- FLUJO MULTIFASICO A TRAVES DE ESTRANGULADORES , 

8.1. - Introducción. 

Se anticipa que la predicción del comportamiento del flujo de mezclas - 
gas— liquido en orificios no es un problema que pueda considerarse re- 
suelto. Existen numerosos estudios sobre este tema y se han encontrado 
varias correlaciones de relacionan el gasto a través del orificio, la 
presión y temperatura antes del orificio y el área de estrangulamiento, 
cuando el flujo es crítico. Algunas de las correlaciones obtenidas es- 
tán basadas en trabajos experimentales y se ajustan razonablemente a los 
rangos probados; sin embargo se desconoce su precisión fuera de esos lí- 
mites. 

En el desarrollo de sus correlaciones los autores han supuesto diversas 
relaciones de presión crítica. Establecer un valor fijo para dicha rela- 
ción implica una simplificación que indudablemente se reflejará en la 
exactitud de las prediciones que se obtengan al aplicar las correlacio- 
nes citadas. Por lo anterior es recomendable que al desarrollar una co- 
rrelación se investiguen las fronteras del flujo crítico y además que las 
relaciones se cumplan para los casos extremos, o Sea : flujo solo de gas o 
solo de líquido. * 

(1-3) 

8.2. - Correlaciones de Gilbert, Ros, Baxendell, Achong. 

A partir de datos de producción, Gilbert desarrollo una expresión apli- 
cable al flujo simultáneo gas-líquido a través de estrangulado res. En su 
trabajo describe en forma detallada el papel del estrangulador en un po 
zo y analiza cual es el efecto sobre la producción de cambios bruscos 
en el diámetro del orificio. 

Tomando como base la relación entre las presiones antes y después de un 
orificio para flujo sónico de una fase, Gilbert recomendé para tener flu 
jo sonico, una relación de 0.588 o menor, entre la presión promedio en 
el sistema de recolección (después del estrangulador) y la presión en la 
boca del pozo (antes de estrangulador). 

Utilizando datos adicionales Baxendell actualizó la ecuación de Gilbert, 
modificando los coeficientes. 

Ros orientó su trabajo al flujo de mezclas con alta relación gas-aceite, 
eií las que el gas fue la fase continua. En su desarrollo llegó a una ex- 
presión similar a Gilbert; pero con coeficientes diferentes. Aparente- 
mente su expresión la comprobó con datos de campo. 

Achong también revisó la ecuación de Gilbert y estableció una expresión 
que valido comparándola con más de 100 pruebas de campo. 

las ecuaciones desarrolladas por los investigadores 
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La forma general de 
citados es : 

Aq r R® 
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Presión corriente arriba, (lb/pg ) 

Producción de líquido, (bl/díai 
Relación gas libre líquido, Cpies3/bl} 

Diámetro del estrangulador , (64 avos de pg) 

Constantes que dependen de la correlación y que toman los v<a 
lores siguientes: 


Correlación 

A 

B 

C 

Gilbert 

10.0 

0.546 

1.89 

Ros 

17.40 

0 . 500 

2,00 

Baxendell 

9.56 

0.546 

1.93 

Achong 

3.82 

0.650 

3 , 8 6 . 

Correlación de Poettmann 

y Beck 

CP y B) (4) 



El modelo de P y B fue establecido a partir del trabajo presentado por 
Ros. La precisión de los resultados obtenidos se comprobó comparándo- 
los con 108 datos medidos. El método fue establecido a partir de un ana 
lisis teórico del flujo simultáneo gas-líquido a velocidad sónica a tra 
ves de orificios y una correlación para el comportamiento p V T de los 
fluidos. No se consideró producción de agua. 

Para que exista flujo crítico se supuso que la presión corriente abajo 
debe ser al menos de 0.55 de la presión en la boca del pozo. Bajo es- 
tas condiciones el gasto en el estrangulador es sólo función de la pre- 
sión corriente arriba y de la relación gas-aceite a condiciones de flu- 
jo. 

La ecuación de P y B es: 

_ 1.549 de 2 / 92 73.6 Pj 

° 73.856 p + Y R \l V, (1 + 0.5 m) 

o g V 1 


) .4513 /r + 0.766 
r + 0.5663 


( 8 . 2 ) 



R. - 

d .- 
c 


A,B ,C. - 


Donde : 

0.00504 ?! Z L (R - R s ) 
r - ~ g 

1 o 


(8.3) 


1 + r 


( P g / '») 


(8.4) 
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Siendo : 

r » Relación gas libre aceite a condiciones de flujo 
V = Volumen especifico del liquido (pies 3 de liq, / Ib de mezcla) 
m = Masa de líquido por unidad de masa de mezcla 

8.4. ~ Ecuación de Ashford^ 

A partir de ion balance de energía y considerando que el fluido se expande 
politropicamente al pasar por el estrangulador , Ashford derivo' una ecua- 
ción que describe el flujo multifásico, bajo condiciones sónicas, a tra- 
vés de un orificio. 

Para compensar la ecuación por las suposiciones incluidas en su desarro- 
llo se introdujo en ella un coeficiente de descarga. Sin embargo al eva- 
luarla, comparando sus resultados con datos medidos en 14 pozos, se en- 
contró que el coeficiente de descarga resultaba muy cercano a la unidad. 

En su derivación Ashford supuso una relación de calores específicos 
k = 1.04 y una relación de presiones, para obtener flujo sónico en el - 
Orificio, de 0.544. 

La ecuación propuesta por Ashford es : 


1.53 d 2 P • 
c 1 

J^(T X + 460) Z ± (R-R s ) + 151 p 

1 ( y + 0.000217 y R +WOR y ) 

lj r ro r g s K w 

0.5 

(B q > W0R)° 

.5 

(T 1 +460)Z 1 (R-R s )+ 111 pj 

( y + 0.00Ó217 y R + WOR y ) 
ro g s r w' 
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Donde : 

d^ = Diámetro del estrangulador (64 avos de pg) 

8.5.- Modelo de Ashford y Piercei — 

Ashford y Pierce establecieron una ecuación aplicable a la región de flu 
jo subcrítico. La validez del modelo se verificó con pruebas de campo - 
diseñadas específicamente. De estas pruebas se obtuvieron valores del - 
coeficiente de descarga en función del diámetro del orificio. Dichos re- 
sultados son: 


Diámetro del 

Coeficiente 

orificio (pg) 

descarga. 

14/64 

1.1510 

16/64 

1.0564 

20/64 

0.9760 
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La ecuación obtenida por los autores es; 


q > 8071,5 C d <x i 
o c ' 


(B q + WOR)‘ 


* k Z 1 
k 


0.00504 C- k ~ _ - x HT r + 460) Z x (R - R g ) (1 - X ■ ) + p '(*> X ) 


( 62 .4 y ro + 0.01353 ^ R +67 WOR ) “ 

(62 * 4 V ro + °* 01353 Y* g Rs + 67 WOR ) 


1 + 0.00504 


(T, + 460) Z, 


(R ~ R c ) X 


k = C P 7 c. 


x - p 2 / p i 


8.6.- Correlación de Omaña R. ' ' 

Omaña desarrollo una correlación Cpara flujo crítico) entre el gasto, la pre 
sión corriente arriba del estrangulador f la relación gas líquido, la 
densidad de los líquidos y el tamaño del orificio. Dicha correlación se 
obtuvo a partir de datos experimentales. En vista de que estos datos estu 
vieron dentro de rangos muy limitados, su aplicación sólo se recomienda ~~ 
para orificios hasta de 14/64 de pg. y gastos máximos de 800 bl/día. 

Las condiciones de flujo crítico se fijaron para una relación de presio- 
nes igual o menor dé 0.546 y una relación gas-líquido mayor de 1.0. 

La ecuación establecida, mediante un análisis de regresión múltiple, es: 
n q 

% = — • ri s“ (8 - 10) 

1.84 ( p / a ) 

i-i * 


n q = 0.263 N -349 Ctr) 3-19 CQ1°‘ 657 (Hj) 1,8 


t> / P T 

g L 
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(.8.13) 


(8.14) 


(8.15) 


La secuencia de cálculo para aplicar la correlación de Omaña puede sinte- 
tizarse en los pasos siguientes: 

1. - Calcular p , p y a a la presión y temperatura existentes antes 

del estrangufador. 

2. - Evaluar N , N , Q y , a las condiciones prevalecientes corriente 

arriba del estrangulador . 

3. - Obtener N con la Ec. (8.11) y q_ con la C8.10) 

q L 


Antes de usar la ecuación (8.10) es conveniente comprobar su validez y 
ajustarla para las condiciones de flujo observadas en ¡un campo. El ajus_ 
te se efectúa introduciendo una constante o coeficiente de descarga del 
orificio (C) , definido por : 

C = gasto medido /gasto calculado. 


8.5 


Flujo Mu£t¿{át>¿co a. tnavte de. EAtMngulcidosieA 


REFERENCIAS . CAP . 8 

1. - Gilbert, W.E. : "Flowing and Gas Lift Well Performance". Drill. 

and Prod. Pract. 1954. 

2. - Ros, N.C. : "An Analysis of Critica! Simultaneous Gas/Liquid 

Flow Through a Restriction and its Application to Flowmetering" 
Appl. Sciences Res., Sec. A. 1960. 

3. - Beggs, H.D. y Brill, J.P. : "Two-Phase Flow in Pipes". Tulsa 

üniversity, 1975. 

4*” Poettmann, F.H. y Beck, R.L.: "New Charts Developed to Predict 
Gas Liquid Flow Through Chokes". World Oil, marzo 1963. 

5. - Ashford, F.E. : "An Evaluation of Critical Multiphase Flow Per- 

formance Through Wellhead Chokes". J.P.T. Ago. 1974. 

6. - Ashford, F.E. y Pierce, P.E. : "The Determination of Multiphase 

'Pressure Drops and Flow Capacities in Down-Hole Safety Valves". 
SPE 5161. AIME 49 th Fall Meeting. Oct. 1974. 

7. - Omaña R. : "Multiphase Flow Through Chokes". SPE 2682. 44th Fall 

Meeting. Sep-Oct. 1969. 


8.6 


9.- COMFORT AMI ENTO DE POZOS FLUYENTES 

9.1. - Introducción. 

Para analizar el comportamiento de un pozo fluyente es necesario consi- 
derar, en forma integral, el sistema de flujo constituido por los si- 
guientes elementos: a) El yacimiento; b) La tubería de producción ; 
c) el estrangulador; d) La línea de descarga. (Fig. 9.1) 

El análisis del sistema de producción se puede efectuar calculando las 
caídas de presión que ocurren en los elementos del sistema, a fin de de^ 
terminar la distribución de presiones en los nodos. (Fig. 9.2) 

El análisis nodal permite determinar la capacidad de producción de un 
pozo; y el efecto del cambio de la T.P., de la línea de descarga, o del 
estrangulador, sobre el gasto. 

La distribución de presiones en el sistema depende del ritmo de produc- 
ción; exceptuando la presión estática del pozo (P ws ) Y la presión de se- 
paración (p ) . . 

-•f 

9.2. - Predicción del comportamiento de un pozo fluyente. 

El método que se verá mas adelante permite obtener las curvas de compor 
tamiento del sistema constituido por el yacimiento, la tubería vertical, 
la tubería horizontal y el estrangulador. Estas curvas (Fig. 9.3) repre- 
sentan la variación de la presión en los nodos y la caída de presión en 
los elementos del sistema, con el gasto. 

Se observa en esta figura que el gasto máximo que puede obtenerse corres^ 
ponde al del flujo sin estrangulador (punto B) . Se aprecia también que 
al estrangular el pozo el gasto disminuye; sin embargo la presión en 
la boca decrece (en vez de continuar aumentando) a partir del punto A, 
en el que se obtiene el gasto óptimo. También se nota el incremento con 
siderable en la caída de presión por la T.P., provocado, como se indi- 
có, por el colgamiento del líquido. 

Para calcular el comportamiento descrito de un pozo, se procede en for- 
ma similar a la indicada para pozos productores de gas. Los pasos son: 

a) Suponer un gasto de aceite. 


b) Con la presión estática y el gasto supuesto, obtener la presión de 
fondo fluyendo. Para hacer esta determinación es necesario usar una 
ecuación que represente el comportamiento del flujo en el yacimiento. 
Esta ecuación puede ser la del índice de productividad, la de Darcy, 
la de Vogel,^ 1 ' u otra. Por ejemplo, se pueden usar las ecua- 
ciones siguientes: 


q 


L 


q o + % 


(9.1) 
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Fig 9.1 Sistema de flujo simplificado 



Fig 9.2 Pérdidas de presión en los elementos 
del sistema de flujo 
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J (P 


P wf ) (l + WOR) 


( 9 . 2 ) 


El índice de productividad a cualquier presión puede obtenerse en, fun- 
ción del índice de productividad inicial (p. s^p ) con la siguiente 
ecuación : 1 


J = J. 


I a oi B oi 


ti B 

^o o 


(9.3) 


Sustituyendo (9.3) en (9.2) y despejando p 


wf * 


-wf 


J. 

i 


( 9 . 4 ) 


(1 + WOR) 


Es evidente que para resolver la Ec. (9.4) se requiere de la predicción 
del comportamiento del yacimiento, ya que los valores de K ro y de WOR - 
son función de la saturación de aceite y de la saturación de agua. 

c) A partir de la presión de fondo fluyendo obtenida se calcula, para el 
gasto supuesto, la presión en la boca del pozo. Este cálculo- se rea- 
liza aplicando el método seleccionado para determinar las perdidas de 
presión en la T.P. El valor de la presión obtenido (p til ) corresponde 

' * - al flujo corriente arriba del estrangulador . Para efectuar este cálcu 
lo, es necesario estimar previamente la relación gas producido-aceite. 

d) A continuación la secuencia de cálculo se reanuda a partir de la pr<s 
sión de separación, para obtener la presión en la boca del pozo co- 
rriente abajo del estrangulador (p- ) * necesaria para transportar el 
gasto supuesto a través de la línea de descarga. 

e) Repetir el procedimiento suponiendo diferentes gastos. 

f) Graficar los valores de las presiones obtenidas (p^, p t ^, p j contra 

los gastos, como se indica en la Fig. 9.3. e 

En dicha figura se observa que cuando p-^h - P e se tiene el gasto máxi- 
mo, correspondiente al flujo sin estrangulador. Los gastos inferiores se 
pueden obtener usando estran guiadores en el cabezal del pozo. El tamaño 
del estrangulador se puede calcular mediante las ecuaciones presentadas 
en el capítulo anterior. 

En la Fig. 9.3 se aprecia también que al ir reduciendo el diámetro de 
los estrangulado res , disminuye el gasto y aumenta la presión en la boca 
del pozo (Pth ) • hasta alcanzar un valor máximo, indicado por el punto 
03 .. Estrangulamientos adicionales provocarán una reducción en la pre- 
sión en la boca, al aumentar las perdidas de presión por la tubería de 
producción. 


Es evidente que la elaboración de figuras como la anterior permite anti_ 
cipar el efecto del cambio de un estrangulador sobre el gasto y la pre- 
sión en la superficie. El manejo inapropiado de los estrangul adores pue^ 


9.3 


Comportamiento de. Pozo ¿ Píayzyvtzs . 



FIG. 9,3.- DISTRIBUCION DE PRESIONES EN UN SISTEMA DE FLUJO. 


de ocasionar la "muerte del pozo". Por ejemplo, un estrangulamiento adi 
cional al obtenido con el orificio correspondiente a 0 3 (cambio dé 0 3 a 
02 ) • originaría un incremento en el colgamiento del líquido y éste, el 
aumento en la carga hidrostática, lo que provocaría una reducción adi- 
cional en la velocidad del flujo. El resultado de esta secuencia, como 
ya se indico, es la precipitación del flujo en un estado inestable (ca- 
beceo) que generalmente conduce a la suspensión de dicho flujo. 

Al repetir el procedimiento de calculo expuesto, considerando Valores 
decrecientes de la presión estática del pozo, se obtienen las relacio- 
nes existentes entre estas presiones y los gastos máximos correspon- 
dientes . (Fig. 9.4) 
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Fig. 9.4.- GASTO MAXIMO OBTENIBLE EN FUNCION DE LA PRESION ESTATICA DEL 
POZO. 

9.3. - Terminación del flujo natural. 

Una aplicación adicional de los procedimientos de cálculo sobre flujo - 
multifásico vertical, es la determinación de la presión estática a la 
cual el pozo dejará de fluir. El procedimiento consiste en graficar 
los valores de la p w f obtenidos a partir del comportamiento del flujo 
en el yacimiento y del flujo por la T.P. (Fig. 9.5) El comportamiento 
de afluencia que se muestra corresponde a las presiones estáticas de 
1200 y 1300 lb/pg 2 . El pozo tiéne una T.P. de 3 1/2 pg. Con una pre- 
sión en la boca de 100 lb/pg 2 , el pozo no fluirá a una presión estáti- 
ca menor de 1250 lb/pg 2 . A una p ws 1150 lb/pg 2 el pozo estará muerto. 
Se advierte que el gasto es de 100 bl Q /dia cuando el pozo deja de produ 
¿ir. Esta situación puede presentarse de un día para otro. En la misma 
figura se observa que con una T.P. de menor diámetro (1.9 pg) el flujo 
natural continuaría por mayor tiempo, hasta que la presión estática se 
abatiera a 900 lb/pg^. 

(4) 

9.4. - Diseño de tuberías de producción y líneas de descarga. 

El procedimiento enunciado anteriormente , permite analizar el efecto del 
cambio de las tuberías de producción y de descarga sobre el gasto. La - 
selección de las tuberías debe basarse en un análisis económico , en el 
que se comparen los incrementos en la producción, al instalar tuberías 
de mayor o diferente diámetro , con la ^inversión adicional que es necesa- 
rio realizar. 
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Fig 9.5 Determinación de la presión esta'tica a la que el pozo 
deja de fluir f3) 
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De este modo pueden determinarse, para cada etapa de la vida fluyente de 
un pozo, cuales son las tuberías necesarias para su explotación óptima. 

Al analizar el efecto del cambio de las tuberías, sobre el gasto máximo 
obtenible, generalmente se obtienen resultados como los mostrados en las 
Figs. 9.6 y 9.7. 

La Fig. 9.6 muestra la variación del gasto máximo al usar líneas de des- 
carga de diferentes diámetros. Se observa que para una tubería de produc 
cion dada, existe un diámetro de línea de descarga para el cual se ob- ~ 
tiene el máximo gasto. Incrementos adicionales en el diámetro de la lí- 
nea de descarga ya no proporcionan mayor producción. : ^ 



FIG. 9.6.- RELACION ENTRE EL GASTO MAXIMO Y EL DIAMETRO DE LA LINEA DE 

DESCARGA, PARA UNA T.P. Y UNA p DADAS. 

■ ws 
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La Fig. 9.7, muestra la variación del gasto máximo al utilizar tuberías 
de producción de diferentes diámetros. Se aprecia que el gasto aumenta 
hasta alcanzar un valor máximo y posteriormente disminuye. 



FIG. 9.7.- RELACION ENTRE EL GASTO MAXIMO Y LA T.P. 

La c omb inación mas adecuada de tubería, se obtiene al analizar diferen- 
tes alternativas y determinar la que permita prolongar al máximo la eta 
pa fluyente del pozo. 

En relación a los procedimientos descritos para determinar la termina- 
ción del flujo natural y diseñar las tuberías de producción, es nece- 
sario indicar lo siguiente: 

1.- Al aplicar cualquier método de flujo multifásico (como los indica- 
dos en las referencias 1 a 9 del capítulo anterior) se obtiene un 
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- gasto optimp de aceite, a partir del cual la presión de fondo aumen 
ta al disminuir el gasto. 

2. - Todos los métodos indican que los gastos óptimos decrecen al dis- 

minuir el diámetro de la tubería de producción. 

3. - Los valores de los gastos óptimos son diferentes para cada corre- 

lación. 


4.- Solo el método de Orkiszewski muestra que la presión de fondo - 
correspondiente a! gasto óptimo, decrece al disminuir el diáme- 
tro de la T.P. Los otros métodos indican la tendencia opuesta. 

9.5.- Optimización de un sistema de producción. 

Los procedimientos de cálculo descritos, aplicados al análisis de un 
sistema de producción dado, permiten identificar los elementos que - 
limitan la capacidad de flujo del sistema. Las principales restric — 
ciones son : 

a) Alta presión de separación en la primera etapa. 

b) Válvulas y conexiones inapropiadas (muy chicas, a 90°, en exceso, 
etc. ) 

c) Tubería de producción inadecuada. 

d) Línea de descarga demasiado larga o de diámetro pequeño, 
ej Línea de descarga común a varios pozos. 

f) Línea de descarga en terreno montañoso. 

g) Daño a la formación en la vecindad del pozo* 

En consecuencia el análisis nodal se realiza para determinar el efec- 
to que, sobre el gasto, tendrían las modificaciones siguientes: 

a) Disminuir la presión de separación . 

b) Eliminar o cambiar válvulas o conexiones inapropiadas. 

c) Colocar separadores a boca de pozo. En este caso se pueden anali^ 
zar dos opciones : , 

1 . - Separar a la presión necesaria para transportar el aceite has ta- 

la central de recolección. 

2. - Separar a baja presión (10-30 lb/pg2) , y bombear el aceite (in 

crementando su presión) hasta la central de recolección. 

*La identificación, prevención y remoción de este tipo de restricción 
al flujo se trata en otras asignaturas. 
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d) Cambiar la tubería de producción. 

e) Cambiar la línea de descarga o instalar una línea adicional. 

f) Instalar un sistema artificial de producción. 

Es evidente la conveniencia de prever desde el principio las restric- 
ciones al flujo, para diseñar el sistema en forma apropiada. Tambien- 
es obvio que la selección de las modificaciones a un sistema y el or- 
den de su aplicación debe basarse en un análisis económico, en el que 
se comparen los incrementos en la producción, al efectuar algún cam* — 
bio, con la inversión adicional que es necesario realizar. 
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lo.- CALCULO DE LA DISTRIBUCION DE LA TEMPERATURA EN TUBERIAS. 

10.1. - Introducción . 

En este capítulo se presenta un conjunto de ecuaciones, para calcular el 
perfil de la temperatura en una tubería que transporta agua, aceite, gas 
o una mezcla de dichos fluidos . v 

La temperatura es un parámetro importante, ya que las propiedades de los 
fluidos dependen de ella y, por consiguiente, las perdidas de presión en 
las tuberías. Generalmente la variación de la temperatura en las tube- 
rías se supone; sin embargo en muchos casos es conveniente calcularla, 
a fin de asegurar resultados mas precisos. Algunos ejemplos son; a) Di- 
seño de oleoductos que transportan aceite viscoso; b) Diseño de tuberías 
submarinas; c) Cálculo de caídas de presión en pozos productores de acei_ 
te volátil o de gas y condensado ; d) Calculó del cambio en la longitud 
de una tubería de producción y en sus esfuerzos, al realizar un trata- 
miento de estimulación, o al producir el pozo .con gastos altos. 

10.2. Perfil de temperatura estable en una tubería con flujo horizontal 
monofásico. 

Sea un tubo de longitud L, del que se considera una sección dx, en la 
que ocurre un cambio de temperatura dT. (Fig. 10.1) 



FIG. 10.1.- ILUSTRACION DE LA TRANSFERENCIA DE CALOR DE 


UN FLUIDO. 
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Al llegar el fluido al punto x, tiene una cantidad de calor Q; pero al - 
pasar a la posición x + dx, pierde cierta cantidad de calor dQ, a través 
de las paredes del tubo. Esta puede ser expresada de la manera siguiente 


dQ 


n d ü dx ( T - Ta) 
12 


( 10 , 1 ) 


En donde:.. 

Q.- Flujo de calor, Btu / dia 
d.- Diámetro de la tubería, pg 

ü.- Coeficiente de transmisión a calor, Btu/dia-pie °F 
T.- Temperatura del fluido, °F 

T a .- Temperatura del medio ambiente que rodea a la tubería, °F 
X.- Longitud, pies 

También puede expresarse dQ de la manera siguiente : í 


dQ - - W f C f dT 


(10.2) 


Que es la pérdida dé calor del fluido al pasar de (x) a (x + dx) ; donde: 

W.- Gasto más ico , lbm/di a 
C^.- Calor específico del fluido, Btu/lb-°F 

Igualando la Ec. (10.2) con la (10.1) 


- WfCf dT 


n d U dx (T - T a )/12 


Agrupando : 


dT 


T - T a 


- a dx 


(10.3) 


Donde : 


n d ü / 12 W f C f 


a = 0.2618 dü/ W f C f 


(10.4) 

(10.5) 


Integrando : 


In 


T - T 
* a 


T_ - T 
1 a 


( 10 . 6 ) 
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Despejando Tx, queda; 

Tx = Ta + (^ - Ta) exp t-ax) (10.7) 

Esta ecuación permite calcular la temperatura estable en un oleoducto. El 
problema para aplicarla es la determinación del valor del coeficiente de 
transmisión de calor. Para oleoductos enterrado^ se ha encontrado experi- 
mentalmente que el valor de U varía entre 4.8 y 14. El valor de C f para el 
aceite varía entre 0.35 y 0.60. 

Para calcular la distribución de la temperatura en un gasoducto, hay que 
considerar el efecto del cambio en la temperatura, resultante por la ex- 
pansión de el gas, al abatirse su presión a lo largo de la tubería. Este 
efecto se conoce como efecto Joule - Thompson. 

La ecuación para un gasoducto es : ^ 

=Jt x - (T a + (Va) (dp/dx) )j exp (-ax) ' 

+ ( T a + ( n/ a) (dp/dx) ) i (10.8) 

Donde : v ' •'"* ;í 

^.-Coeficiente de Joule - Thompson, definido por; (dT / dpi. Su valor 
puede obtenerse de la Fig. 10.2, derivada de la referencia J2. ' 

dp/dx.- Gradiente de presión, lb/pg 2 /pie 

Para resolver la Ec. (10.8) es necesario previamente calcular el valor de 
dp/dx, según se vio en el capítulo 4. f t i 

Ejemplo.- Se desea calcular la temperatura del gas, en un gasóducto ente- 
rrado, del que se tiene la siguiente información; 

x = 30 millas = 158310 pies V 

d = 12 pg 

Cf g ■ 0.8 Btu/lb - °F 

q » 45 x 10 6 pies 3 a c.s./dia 

Vg = 0.698 

T^ = 150°F (Temperatura inicial, a x = 0) 

T a = 70°F 

ü = 9.6 Btu/dia pie 2 °F 

dp/dx = - 8.5276 x 10”^ lb/pg 2 /pie . 

0.1 °F/lb/pg 2 
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Y„ ^ 3 ^P^ - e . ) g 0.0764 (lb/pie 3 ) a 

y (Ib /pie* 3 ) a 

W = (0.0764) (45 x 10“) (0.698) 

W * 2.4 x 10 6 Ib /dia 

m 

b) Cálculo de a; según la Ec. (10^5): 

»■ 

a = (0.2618) (12) (9.6) / (2.4 X 10 6 ) (0.8) 
a = 1.5708 x 10 -5 . , 

c) Otros factores de la Ec (10.8) son: 

exp ( -ax) = exp ({- 1.S708 x 10 -5 ) (158310)) = 0.08318 

í'I/a) (dp/dx) = (0.1/1.5708 x lo" 5 ) (-8.5276 x lo' 4 ) = - 5.4288 


Solución. - 
a) Cálculo de W. - 


W ( Ib ] 

\diá )~ 


q g (: 


pies a a c. 


dia 


W = 0.0764 q v 

g y g 


d) Sustituyendo en la Ec. (10.8): 

: T- = (150 - (70 - 5.4288)) 0.08318 + (70 - 5.4288)) 


71. 68°F 


La solución rigurosa de la Ec. (10.8) comprende la determinación previa 
de t), de acuerdo con : f?) 


8 T 
8P 


SE 

~rr 

SE 

ST 


(10.9) 


Donde H es la entalpia específica, cuyos valores se calculan mediante la 
aplicación de una ecuación de estado, a diferentes presiones y tempera- 
turas. Las derivadas de la entalpia con la presión y con la temperatura 
se obtienen del cálculo anterior. Como H y en consecuencia rj dependen de 
la composición del fluido , la presión y la temperatura, para resolver 
la Ec. (10.8) es necesario utilizar un proceso iterativo, ya que 
dp/dx depende tairbién de la temperatura. El proceso iterativo consiste en 
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a) .- A partir de un punto de presión y temperatura conocidos (P,, T ) , 

considerar un incremento de longitud Ax y suponer una caída ie - 
presión y una caída de temperatura. 

b) .- Determinar la presión media y la temperatura media del tramo de tu- 

bería considerado, y calcular las propiedades del fluido a: p ,T ; 
P 2 r T 2 y p, T. 11 

c) .- Calcular p 2 con la Ec. (4.16) y dp/dx 

d) .- Obtener rj con la Ec. (10.9) y T con la (10.8) 

e) .- Repetir el procedimiento hasta que los valores de p y T coincidan 

con los supuestos. 22 

10.3.- Perfil de temperaturas en flujo horizontal multifásico. 

La siguiente ecuación puede emplearse para calcular la temperatura cuan- 
do se tiene flujo multifásico: (2) 

+ (T^ - T ) exp (-ax) 

a=0.2618üd/(W C + W C. + Vi C, ) (10.10) 

g fg o fo w fw 

En esta ecuación Cfg, C fo y C fw , son los calores específicos del gas, - 
del aceite y del agua. Los gastos másicos pueden obtenerse con: 


II 

p g % 

CR 

- «s 

) Bg 

(10.11) 

w o ■ 

5.615 

P o 

V 

B 

O 

(10.12) 

W w = 

5.615 

P w 


B 

w 

(10.13) 


La dificultad obvia en resolver éstas ecuaciones ha conducido al uso de 
valores de Wg, W Q y W w a condiciones superficiales, asi como al empleo de 
valores típicos de Cfg = 0.8, 0f o * 0.4 y de = 1.0 Btu /lb m °F 

El valor de Cfg puede obtenerse de la Fig. 10.3 

10.4.- Perfil de temperaturas en tuberías de producción. 

La temperatura a lo largo del pozo puede determinarse, en forma senci- 
lla y precisa, mediante el procedimiento propuesto por Romero Juárez.' * 

> ^ (4) 

En el calculo se emplea la ecuación de Ramey: 

T (Z, t) = g Q Z + T es - g e A + (T fs + g Q A - T es ) exp C-Z/A) CIO. 14) 
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Calor específico (BTU/lb-°F) Cfg 



Fig 10.3 Calor epectfico del gas natural con densidad relativa de 
0.6 a 0.7 
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8 ” r ti ü 7 (10.21) 

La nomenclatura usada en estas ecuaciones es : 

A = función definida por la Ec. (10.15), pies 

B = cantidad definida por la Ec. (10.21) 

= calor específico del fluido, Btu/lbjQ - °F 
<* ce * diámetro exterior de la T.R. , pg ' ^ Vh- 

d ci = diámetro interno de la T.R. , pg 

d fce = diámetro externo de la T.P. , pg 

dti = diámetro interno de la T.P,, pg . 

F = factor definido por la Ec. (10.20) 

f(t)= funcipn del tiempo de la conducción de calor 
g e = gradiente geotérmico, °F/pie 

li = coeficiente de transferencia de calor de la película interior, 
Btu/dia - pie 2 - °F 
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“han 

*Sie 

*h. 


bh 


fs 
ü , 


= gasto de inyección, bl/dia 

_ conductividad térmica en el espació anular , Btu/dia- pie - °F 
= conductividad térmica de la tierra, Btu/dia - pie - °F 
* conductividad térmica del acero, Btu/dia - pie - °F 
= radio interior de la T.P., pies ¡, 

= tiempo , dias 

= temperatura en el fondo del pozo, °F 
» temperatura del terreno en la superficie, °F 
= temperatura del fluido en la superficie, °F 

2 

= coeficiente de transferencia de calor total, Btu/di.a-pie -°F 
= gasto másico, lb m /dia 
= amplitud del espacio anular , pies 
= espesor de la T.R. , pies 
= espesor de la T.P., pies 


q = densidad del fluido lb m /bl 


10. 4.1.- Temperaturas en un pozo inyector de agua. 

Usando valores típicos de - 350 lb/bl, - 1 Bt u /lb m - °F , 

K he = 33.6 Btü/dia-pie - °F, K hs = 600 Btu/dia - pie - °F y 

K han ~ 9-432 Btu/dia-pie- °F (suponiendo el espacio anular lleno dé agua) 

las ecuaciones (10.2.0) y (10.21) se convierten en: 

( 10 . 22 ) 

(10.23) 

EÍ valor de U, despreciando la resistencia al flujo de calor de los dos 
primeros términos de la Ec. (10.16), ya que son pequeños comparados con 
la resistencia al flujo de calor del espacio anular, queda: 

U - 226.4 / (d ci - d te ) (10.24) 


F = 1.658 (1 + B f (t) )/B 
B = d fci U / 806.4 


Sustituyendo en la Ec. (10.21) 

B - 0.281 dti / (d ci - d te ) 


(10.25) 
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Para pozos inyectores a través de la T.R, (sin T.P. ) , d te = d c ¿, por lo 
que : 

A » 1.658 i f(t) (10.26) 

Cuando el espacio anular contiene aceite con: 

K¡ lf = 1.896 , V ro - 0.8 y = Q.3 c.p. , el valor de B, es: 

B * 0.0564 d ti /(d ci - d te ) (10.27) 

Ejemplo: Se desea calcular la temperatura de un pozo inyector de agua, 

a 14 000 pies, para un gasto de 2 400 bl/día, después de 50 días de in- 
yección. El gradiente geotérmico es de 0.019 °F/pie. Otros datos son: 
Tfs = 60 ®F , T es « 70°F, d ci = 5.675 pg. , d ce =6.625 pg. 

Solución: ; , 

Y = (552) (50) / (6.625) 2 = 628.8359 

f (t) = 3.5732 

De la Ec. (10.26) 

A = (1.658) (2400) (3.5732) =14218 
Sustituyendo en la Ec. (10.14) 

T (Z, t) = (0.019) (14000) + 70 - (0.019) (14218) 

+;( (60) + (0.019) (14218) - (70) ) exp. (-14000/14218) 

T = 163 °F 


10.4.2.- Temperaturas en un pozo productor de aceite y gas 


(5) 


Los pozos productores con gastos altos pueden sufrir elongaciones con si_ 
derables en sus T.P. También es importante estimar lá temperatura en la 
boca del pozo para definir el procesamiento de los fluidos producidos. 


La Ec. (10.14) , para un pozo productor, se convierte en: 
T (Z, t) * T bh + g e ( A (1 - exp (-Z/A) ) - Z ) 

El valor de A se calcula con la Ec, (10.15), donde: 


= W 

f m . 


(10.28) 


(10.29) 


Para usar la Ec. (10.15) es necesario determinar previamente el valor de 
Cf , que puede obtenerse ponderando los calores específicos de cada fase 
con sus respectivos gastos másicos. 


(W + W C- + W C_ ) / Ví 

v g fg o fo w fw ■ ni 


(10.30) 
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Cuando el espacio anular contiene aceite, el valor de ü puede aproximar- 
con : 

• U = 45.50 / (d ci - d te ) (10.31) 

En esta ecuación se considera que la conductividad térmica del aceite 
es de 1.896 Btu/dia-pie-°F 3 

Ejemplo: Calcular la temperatura en la boca de un pozo que produce 
15 000 blo/dia, al cabo de 30 días; si y ro = 0.835, V = 0.75, 

R = 2000 pies 3 /bl. El yacimiento está a 9031 pies de profundidad y 
T bh es de 212 °F. El gradiente geotérmico es de 0.019 °F/pie. Los diám^ 
tros de las tuberías son: d ce =9.626 pg. d c ¿ = 8.535 pg. d t ¿ = 6.094 pg 
y ^te = 7 -000 pg. El espacio anular contiene aceite. El calor especí- 
fico de los fluidos producidos es de 0.55 Btu/lb -°F. La producción de 
agua es nula. 

Solución: 

Sustituyendo en la Ec. (10.15) , los siguientes términos: 
ü = 45.50/ (d ci - d te ) 


\e= 33 ' 6 


W = q M 
f o 


K =1.896 

nan 

■■■ ■ r . . " 

Se tiene: 

ft= V" C, W el . - V_ (33.6 + 1.896 d t . f (t) / (d c . - á te ) ) (10 32) 
400.27 d 

ti 

El valor de f(t) de acuerdo con las Ecs. (10.17) y (10.18) , para t = 30 
días y d ce = 9.626 pg, resulta: ■ 

f (t) = 3.0025 

El valor de M, según la Ec. (5.16) , es: 

M = 350.5 y + 0.0764 R y 
r ro g 

.*. M = 407.27 Ib / bl a c.s. 

m o 
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Sustituyendo valores en la Ec. (10.3 2) : 

A a (15 000) (407.27) (0.55) ( (33.6)+(1.896) (6.094) (3.0025) /(l. 535) ) 
(400.24) (6.094) / (1.535) 

A = 118834.7 

Sustituyendo en la Ec. (10.28): 

T > 212 + 0.019 (118835 (1 - exp (-9031/118835) ) - 9031) 

T , = 205.6 °F 
tn 

El mismo ejemplo, para q Q = 1000 blo/dia, 
da: 

T th * 142.8 °F 

10.4.3.- Efecto Joule - Thompson en flujo multifásico. 

Este efecto se presenta solo cuando la presión es menor que la de satu- 
ración. Cuando el gas se libera en cantidades apreciables, su efecto de^ 
be ser considerado en la caída de la temperatura. El abatimiento en la 
temperatura, debido a la liberación y expansión del gas, puede calcular 
se con : (5 ) > 

AT b = k’ (Z - Zh) (10.33)" 

Donde: 

k' = constante de abatimiento de la temperatura por el efecto Joule- 
Thompson , °F/pie 

* Longitud a la que se tiene la presión de saturación del aceite, pies 

Para considerar el efecto descrito en los cálculos, basta introducir 
A T b en las ecuaciones (10.10) y (10.28). 

El valor de k' , para pozos productores en Arabia con una relación gas- 
aceite de 540 pies 3 /bl, resultó de 0.0015 °F/pie? sin embargo para al- 
tas relaciones gas aceite puede ser mayor dicho valor. 

10.4.4*- Abatimiento de la temperatura en las tuberías conductoras ver- 
ticales (risers) utilizadas en los sistemas de producción sub- 
marinos. 

El ritmo de transferencia de calor en "las tuberías submarinas es gene- 
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raímente mayor que el existente en tuberías subterráneas. Esto se debe a 
las corrientes de convección del agua, que disipan el calor mas rápida- 
mente que en la conducción para tuberías enterradas o la convección forza- 
da de aire sobre lana tubería superficial. 

Las tuberías conductoras verticales generalmente están descubiertas, - 
mientras que las líneas de recolección submarinas están recubiertas con 
concreto, para compensar el efecto de flotación. El efecto aislante de 
calor, del concreto, debe incluirse en los cálculos del abatimiento en 
la temperatura. 


El perfil de la temperatura se calcula con la Ec. (10.10) , cuando se ti£ 
ne flujo multi fásico. El problema para aplicar esta ecuación es el cálcu 
lo de U, que está dado por: 



(10.34) 


Donde : 


= espesor del recubrimiento de concreto, pies 
X fc = espesor de la tubería, pies 

K^ c = conductividad térmica del recubrimiento de concreto o del mate- 
rial aislante, Btu/dia - pie - °F. 

2 

h = coeficiente térmico de la película exterior, Btu/dia - pie-°F 
o 


Generalmente h puede despreciarse. Cuando el recubrimiento de concreto 
es mayor de 4 pg., también son despreciables K^g y h Q , La conductividad 
térmica del concreto puede variar substancialmente con el contenido de 
humedad y su integridad estructural. La conductividad térmica de la pe- 
lícula exterior es función de 1 a velocidad del fluido ambiental, normal 
a la tubería. 

El valor de h Q puede obtenerse mediante la aplicación de diversas co- 
rrelaciones. Aquí se verá solamente la establecida por Fand, que es:^ 7 ' 

N u = ( 0‘.35 + 0.56 Re ° -52 ) Pr °' 3 (10.35) 

En donde todas las propiedades se determinan a la media aritmética de 
las temperaturas del flujo libre y de la pared. 

En esta ecuación : 

N u = numero de Nusselt de la película exterior. ^ 

R e = numero de Reynolds de la película exterior. 

P r = numero de Prandtl de la película exterior. 
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Los valores de estos números adimensionales se pueden obtener con las - 
siguientes ecuaciones: 


Nu = 


o d 
12 K_ 


(10.36) 


R = 124 d Vf P f 

e V f 


(10.37) 


P = 58.06 
r 


c f H 


(10.38) 


Ejemplo: Calcule el valor de U para una tubería conductora vertical des 

cubierta, de 10 pg de diámetro exterior, sobre la cual se tiene una co- 
rriente de agua con una velocidad dé un pie/seg. La temperatura del agua 
es.de 50°F . El espesor de la tubería es de 0.02604 pies. La temperatura 
del fluido transportado en el interior de la tubería es de 190 °F. 

Solución: 


La temperatura media del fluido es: 

+ 

2 


r . i 2g_i j j o; = 120 . F 


Las propiedades del agua a esta temperatura , son : 


u = 0.65 c.p. 
r w * 

K = 8.4 B tu/di a - pie - °F 


C f = 1 Btu/lbm - °F 


P- = 61.5 Ib 

, 3 
pie 


El número de Reynolds es : 


R . JimMmi- - 117323 
e 0.65 


El número de Prandtl es: 

(58.06) (1) (0.65) 


P = 

r 


(8.4) 

Sustituyendo valores en la Ec. (10.35) 


4.4927 


N 


h n d 


u 12 ^ 

( 380.756) (12) (8.4) 
( 10 ) 


( (0.35) + (.0.56) (117323) °* 52 ) (4.4927) 


0.3 


= 3838 
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Sustituyendo en la Ec. (10.34) , en la que se considera un valor de K. 
600 , se tiene : * ls 


k hs 


Sustituyendo : 
1 


3838 


0.02604 

600 


ü 35 3290 Btu/dia - pie - °F 
10.5.- Coeficiente de transferencia de calor. 

Según se indico, el cálculo del coeficiente de transferencia de calor 
es en sí un problema importante en la determinación del perfil de la 
temperatura en una tubería. Debido a la variedad de materiales que ro 
deán a las tuberías, puede existir una mezcla compleja de perdidas de 
calor. Por ejemplo, el espacio anular de un pozo revestido puede con- 
tener cemento, gas, agua, lodo o aceite. El valor del coeficiente de 
transferencia de calor puede modificarse además por la presencia de 
parafina, incrustaciones, recubrimientos anticorrosivos, etc. El ter- 
mino h es extremadamente complejo en flujo multifásico, ya que depen- 
de del patrón de flujo, además de los parámetros normalmente acepta- 
dos en flujo de una fase. 


El valor de h puede estimarse usando la siguiente ecuación estableci 
da por McAdams. ^ 


N u 


hd 

*hf 


0.023 Re°* 8 Pr°* 4 


(10.39) 


Donde los valores de N u , Re y Pr, corresponden al del flujo en el int<e 
rior de la tubería* 


La ecuación anterior, expresada en unidades prácticas, es: 
hd 

= U.Ü23 t 

p. 


12 K = 0.023 (HÍÍILA. i 0 '® ( - 5a -° 6 K Cf )°' 4 (10.40) 

12 K hf ü f \f 


El coeficiente de transferencia de calor cuando se tiene una tubería 
enterrada puede obtenerse con la siguiente ecuación: 


■♦■Jfcr + * 


ln (48Z/d) 


«he *1,3 24 Khe /d 


(10.41) 


Donde : 
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x r ® espesor del recubrimiento , pies 

Z = distancia entre la superficie del terreno y el centro de la 
tubería# pies. 

- conductividad térmica del suelo# Btu/dia-pie °F 
d = diámetro internó de la tubería, pg. 


En la Tabla 10.1 se muestran los valores de Kft e de acuerdo con las ca- 
racterísticas del terreno. 


TABLA 10.1 VALORES DE 




( 10 ) 


Tipo de suelo 

Grado de humedad 

K he 



(Btu/dia pie 

Arenoso 

Seco 

10 - 20 

Arenoso 

Húmedo 

22 - 29 

Arenoso 

Saturado 

53 - 63 

Arcilloso 

Seco 

9 - 14 

Arcilloso 

Húmedo 

19 - 24 

Arcilloso 

Saturado 

29-44 


Otros valores típicos de materiales usados como recubrimiento, son: 

Btu/dia pie °F 


( 11 ) 


Espuma de uretano 

Poliestireno 

Concreto 


0.22 

0.54 

24.00 


10.10 
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10.17 


NOMENCLATURA: 

Las unidades expuestas aquí corresponden a las de las ecuaciones que se 
usan en el texto para la resolución de problemas prácticos. Las ecua- 
ciones básicas están en unidades consistentes. 


d ci 

%. 

d te 

d ti 

D 

E 

E 


tp 

f 

o 

f 

w 

f(t) 


2 

area, pg 

área anular, pg 2 

3 ; 3 

factor de volumen, pie /pie 

coeficiente «de descarga del estrangulados, adimensional. 

2 

compresibilidad, pg /Ib 
calor específico del fluido, Btu/lb °F 
diámetro, interno de la tubería,pg 
diámetro del estrangulador , 64 avos de pg 
diámetro exterior de la T. R. , pg 
diámetro interior de la T. R. , pg 
diámetro hidráulico, pg 
diámetro exterior de la T.P. , pg 
diámetro interior de T.P. , pg 
profundidad, pies 

energía por unidad de masa, lb^-pie/lb^ 

factor de eficiencia, adimensional 

^ 2 

perdidas de presión por aceleración, lb f /pg 
factor de fricción 

factor de fricción del diagrama de Moody para tuberías lisas 
factor de fricción para dos fases 
relación ace i te -líquido a c. esc. 
relación agua-líquido a c. esc. 

función del tiempo en la conducción del calor en regimen variable 

„ 2 

aceleración de la gravedad, pies/seg 


g 


factor de conversión en la 2a. ley de Newton, Ib -oie/lb -seg 2 
gradiente geotérmico, °F/pie 
elevación , pies 

coeficiente térmico de la película interior , Btu/día-pie 2 - °F 
coeficiente térmico de la película exterior, Btu/día-pie 2 - °F 
entalpia específica, Btu/lb 

m 

fracción del volumen de la tubería ocupada por líquido 
("colgamiento”) , 

gasto de inyección, b 1/día 

' ■, ' • ' o 

índice de inyectividad, bl/día/lb f /pg 

2 

índice de productividad , bl/día/lb f /pg 
relación de calores específicos 

conductividad térmica en el espacio anular, Btu/día-pie °F 

conductividad térmica del recubrimiento del material 

conductividad térmica del suelo, Btu/día-pie °F 

conductividad térmica del fluido, Btu/día-pie °F 

conductividad térmica del acero, Btu/día-pie °F 

constante de abatimiento de la temperatura por el efecto Joule- 
Thompson , ° F/pie . 

longitud, millas 

línea de descarga 

longitud equivalente, millas 

logaritmo natural 

logaritmo decimal 

masa. Ib 

m 

peso molecular. Ib /mole -Ib 
m 

masa asociada a un barril de aceite. Ib /bl a c.s 

m o 

número de moles, mole-lb 


gv 


FR 


Reg 


ReL 


pe 


pr 


th 


wf 

P 

ws 

q' 


número de la viscosidad del líquido, adimensional 

numero de Nusselt de la película de fluido, adimensional 

numero de la velocidad del gas, adimensional. 

número de la velocidad del liquido, adimensional 

número de Froude, adimensional 

número de Reynolds, adimensional 

número de Reynolds del gas, adimensional 

número de Reynolds del líquido, adimensional 

presión, lb/pg 2 

2 

presión de burbujeo, lb/pg 

2 

presión corriente abajo del estrangulador , lb/pg 
2 

presión de bombeo , lb/pg 

número de Prandtl de la película de fluido, adimensional 

2 

presión de separación, lb/pg 

2 

presión a las condiciones estándar lb/pg 
> 2 

presión pseudocritica, lb/pg 
presión pseudoreducida 

2 

presión en la boca del pozo, lb/pg 

2 

presión de fondo inyectando, lb/pg 

2 

presión de fondo fluyendo, lb/pg 

2 

presión estática (de fondo cerrado) . lb/pg 

gasto de producción a c. esc. b 1/día 

gasto, bl/día 

gasto óptimo, bl/día 

flujo de calor, Bt u/día 

relación gas libre -aceite a c.e. 

radio interior de la T.P. , pies 


Mi 



R 

R 

e 


sb 

R 

sr 

Rs_ 


T 




T 

pr 


T.R. 


U 

v 

V 

m 



radio de drene, pies 
radio del pozo, pies 

2 3 

constante de los gases reales, (lb/pg -pie ) (°R mole-lb) 
relación gas-aceite, pies 3 /bl 

numero de Reynolds de la película de fluido, adimensional 

radio hidráulico, pg 

relación de solubilidad, pie 3 /bl 

relación de solubilidad corregido, pie 3 /bl 

relación de solubilidad a o , pie 3 /bl 

3 

relación de solubilidad real, pie /bl 

relación de solubilidad obtenida del laboratorio, pie 3 /bl 

relación de solubilidad del gas en el agua, pies 3 /bl 

tiempo, días 

temperatura, °F 

temperatura de separación, °F 
* - 

temperatura a las condiciones estándar, °F 
temperatura en el fondo del pozo, °F 
temperatura de la roca en la superficie, °F 
temperatura del fluido en la superficie, °F 
temperatura pseudocrí tica, °F 
temperatura pseudoreducida, °F 
tubería de producción 
tubería de revestimiento 

2 

coeficiente de conductividad térmica, Btu/día-pie °F 

velocidad, pie/seg 

velocidad de la mezcla, pies/seg 

velocidad real del líquido, pies/seg 


5 * 


SL 


sg 


V 

V 

V 


z 

z 

Ap/AL 

v g 

V 

^gd 

r g£ 

r gs 

Vgp 

*ro 

£ 

n 

w 

0 


velocidad real del gas, pies/seg 
velocidad superficial del líquido, pies/seg 
velocidad superficial del gas, pies/seg 
volumen, pies 3 
volumen específico, pie 3 /lb 

m 

volumen de la mezcla a. c. de esc. por barril de aceite 
producido a.c.s., pies 3 /bl 

gasto másico. Ib /día 
m 

gasto másico. Ib /seg 
m 

relación agua-aceite a c.s. bl /bl 

w o 

amplitud del espacio anular, pies 

espesor de la T.R. , pies 

espesor de la T.P., pies 

factor de compresibilidad del gas 

profundidad bajo el nivel dé la superficie, pies 

^ 2 

gradiente de presión, lb^/pg /pie 
densidad relativa del gas (.aire- 1.00) 
densidad relativa del agua 
densidad del aceite, °ÁPI 

densidad relativa del gas disuelto (aire=1.000) 
densidad relativa del gas libre (aire=1.000) 

2 

densidad relativa del gas a una p^ - 100 lb/pg man. 

densidad relativa del gas a p y T 

s s 

densidad relativa del aceite producido (agua=1.000) 
rugosidad, pg 

2 

coeficiente de Joule-Thompson , °F/lb^/pg 

relación de gastos másicos (w /w ) 

g b 

ángulo de la tubería con la horizontal, grados 


Pg 

Po 

P w 

?<* 

P m 

e n£ 

,1 

> u o 


- w ob 

w om 


L 

r í 


Af 
A 
Ah 
AL ; 


A p 


Av 


A w 


densidad del gas, Ib /pie -3 
m 

densidad del aceite, Ib /oie 3 
m " 

densidad del agua. Ib /oie 3 
tn 

3 

densidad del aceite a o, , Ib /pie 
~b m 

densidad de la mezcla, Ib /pie 3 
m 

densidad de la mezcla sin resbalamiento, Ib /oie 3 

m / 

densidad del líquido. Ib /pie” 3 
m 

colgamiento sin resbalamiento 

viscosidad del aceite, cp 

viscosidad del gas, cp 

viscosidad del aceite a o, , cp 
* b 

viscosidad del aceite muerto, cp 
viscosidad del líquido, cp 
viscosidad de la mezcla, cp 

viscosidad de la mezcla sin resbalamiento, cp 
tensión superficial . del aceite, dmas/cm 

tensión superficial de la mezcla, dinas/cm 

2 

coeficiente de Joyle -Thompson . °F/lb/pg 

densidad del fluido, Ib /bl 
m 

incremento finito 
incremento de elevación, pies 

incrementó de longitud, pies 

„ . . ^ 2 
caída de presión, lb/pg 

incremento de velocidad, pies/seg 

pérdida de energía, Ib -pie/lb 
r " ra 


*j 


SUBINDICES 


a aire 

b burbujeo 

c.esc condiciones de escurrimiento 

es condiciones atmosféricas 

f fricción 


g gas 

L líquido 

m masa 

m mezcla 

^ ns sin resbalamiento 

0 aceite 

s separación 

w agua 

1 corriente arriba del estrangulador 

2 corriente abajo del estrangulador 


vi i 


CONSTANTES Y FACTORES DE CONVERSION 


Condiciones base (c.s.) . 14.7 Ib/pg y 60 °p 

Temperatura absoluta, correspondiente a 0°F . 460° R 

Peso molecular del aire seco ...... .i. 28.97 lb/mole-lb 

Volumen de 1 mole - Ib de gas a c.s.. . ... 379.4 pies 3 ^ 

Densidad del agua a c.s 62.428 lb/'pie^ 

Densidad del aire a c.s. . . . . . . . . . ... . 0.0764 lb/pie' 

Carga hidrostática de 1 pie de agua a c.s. ... . 0.433 lb/pg^ 

R .............. . .... .f, 10.73 ( lb/pg 1 2 -pie 3 ) / ( °R . mole-lb) 


Conversiones 


Longitud 


1 pié = 30 . 48 cm 
1 milla = 5277 pies 


141.5 

131.5 + °API 


y= °api = 


1 atm = 1.033 kg/cm 
1 kg/cm 2 =14.223 lb/pg 2 
1 atm =14.696 lb/pg 2 abs 


i gr /cm 3 = 62.428 Ib /pie 3 
1 gr/Cm 3 = 350.63 lb/bl 
1 Ib /pie 3 = 5.6166 lb/bl 


°F = 1.8 °C + 32 
°C = 5/9 ( °F - 32) 


°R = °F + 460 


1 bl = 42 gal 
1 bl = 5.6146 pie 
1 m 3 = 6.2896 bl 
1 m 3 = 35.314 pie 
1 bl =158.987 lt 


1 Ib = 453.59, gr 

1 kg = 2.2046 Ib 


Transferencia de calor 


1 Btu = 1.055 x 10 joules 
1 watt = 3.4121 Btu/hr 
1 v/att/m = 1.0403 Btu/hr pie 
1 watt/m^ °c = 0.1761 Btu/hr pie 2 


Viscosidad 
1 cp = 6.7197 x 10" 


Ib / pie-seg 


T cp = 0.04031 Ib /pie-min 



